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Summary 


A  method  of  predicting  creep  in  bending  from  data  on  creep   in  tension 
has  been  derived  and  applied  to  creep  of  a  canvas   laminate.     Both  creep 
deflections  and  stress  distribution  were  derived.      Th-e  deflections  compared 
favorably  with  test  data.      It  was   shown  that  the   stress  distribution  re- 
mained constant  during  creep   in  bending  when  creep  in  tension  and  con- 
pression  ¥\rere  equal  and  the  coefficient   of  the  time- dependent  term  "was 
equal  to   the  time- independent  term.      Methods   of  determining  creep  deflec- 
tions  of  beams  having  non-uniform  bending  moments   are  described. 


Conclusions 

The  form  of  creep   equation  employed  seems  adequate  to  describe  the 
tension  creep  of  a  grade-C  canvas  laminate.      The  close  agreenent  between 
creep  in  bending  predicted  by  this    equation  and  actual  data  lends    support 
to  the  possible  usefulness   of  the  equation.     Further  evidence  is   pro- 
vided by  the  agreement  betrvveen  the  stress-strain  curve  derived  from  the 
same  creep  equation  and  actual  test  data. 
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STATUS  OF  THE  PROBLEM 

Tensile  Creep 

numerous  attempts  have  been  made  to  express  creep  behavior  by  a 
mathematical  expression  either  emperically  or  theoretically.  Many  in- 
vestigators have  assumed  the  time  dependence  of  creep  at  constant  load  to 
be  a  linear  function  of  time  (1-4)  and  have  examined  the  rate  of  creep  as 
a  function  of  stress  and  temperature  -  usually  only  stress  has  been  con- 
sidered.  Some  of  the  types  of  equations  XYhich  have  been  employed  to 
describe  the  relationship  betrrireen  creep  rate  v  and  stress  o  are  as  follo'vs: 

V  =  Ba  (viscous)  Eq.  1 

V  -  Bo     (pov/er)  ^q,  2 

Ba 

V  =  A'e   or  In  V  =  A  +  Ba  Eq.  3 

(exponential  or  logarithmic) 

V  =  B  sinh  —  (hyperbolic  sine)        Eq,  4 

o 

where  a.  A,  A',  B,  and  a     are  constants. 

Of  these  equations,  the  hyperbolic  sine,  Eq.  4,  has  been  shoivn  to  have 

fundamental  significance  in  terms  of  molecular  behavior  (5).   It  should 

be  noted  that  Eqs,  1  and  3  are  approxim.ations  of  Eq,  4  if  —  is  either 

o 
small  or  large  respectively. 

Other  investigators  (6-14)  observing  that  creep  was  usually  not  a 

linear  function  of  time,  expressed  the  strain,  e,  versus  time,  t,  relation 

non-linear ly  by  equations,  such  as  the  following: 

e  =  C  +  D  In  t  +  Et  Eq.  5 

See  references  (7,  8) 

e  =  C  +  Dd-e"^"^)  +  Et  Eq.  6 

See  reference  (14) 

e  =  e  +  mt  Eq,  7 

o 

See  references    (6,    9,    11-13) 


+ 

lumbers  in  parentheses  refer  to  the  list  of  references  appended  to  this 
paper. 


2. 

where  b  and  n  are  constants  and  C,   D,   E,    e   ,   and  m  are  functions  of  stress, 

ViTien  the  coefficients   e     and  m  in  Eq.    7  are  expressed  as  hyperbolic 
sine  functions   of   stress,    Eq.    7  becomes    (13,    15): 


■   1-  cy     .      ,,n 
e  =  e     smh  —  +  m't 

O  Q 

S 


smh  — 
o 
m 


Eq.    8 


vjhere  e  ',  m',  a,  a,  and  n  are  constants.  From  Eq,  8  the  creep  rate  v, 
O         £    21 

as  obtained  by  differentiation,  is: 

V  =  m'nt    sinh  —  Eq,  9 

G 

m 
and  is  a  function  of  time  rather  than  a  constant  as  often  assumed. 


Bending  Creep 

The  first  analysis  of  creep  in  bending  was  reported  by  McCullough  (16) 
in  1933,   The  results  of  his  tests  of  lead  beams  virere  in  reasonable  agree- 
ment with  the  deflections  calculated  from  observed  minimum  creep  rates  in 
tension  and  compression.   Linearity  of  strain  with  distance  from  the  neutral 
axis  was  demonstrated,  that  is,  plane  sections  remained  plane  even  during 
creep, 

Tapsell  and  Johnson  (2)  also  demonstrated  for  lead  beams  that  the 
strain  remained  a  linear  function  of  distance  from,  the  neutral  ar.is.   It 
followed  from  this  observation  and  was  sho'vm  by  the  data  that  the  strain- 
time  curves  at  every  position  (fiber)  of  the  beara  have  a  similar  shape. 

Creep  of  beams  of  aluminum  and  polystyrene  plastic  were  investigated 
by  Marin  and  Zwissler  (17)  and  Marin  and  Cuff  (l8).   The  observed  minimum 
creep  rates  in  bending  were  compared  with  bending  creep  rates  computed  from 
tension  creep  tests  assuming  a  linear  dependence  of  creep  on  time  and  a 
power  relation  between  creep  rate  and  stress.   The  theory  did  not  permit 
computation  of  deflection.  Both  the  tension  and  bending  creep  curves 


3. 

depart  considerably  from  the  assumed  linear  dependence  on  time.   In  a 
discussion  Findley  (18)  presented  data  on  creep  of  polystj^rene  -.vhich  in- 
dicated that  the  power  function  of  stress  vra.s  not  satisfactory  for  small 
stresses. 

Creep  of  an  imaginary  copper  beam  v;as  examined  by  Popov  (19)  using 
several  creep  rate  vs.  stress  relations.  liis  examination  indicated  that  the 
stress  distribution  in  a  beam  changed  with  time  and  tended  tov;ard  a  steady- 
state  distribution  of  a  non-linear  shape.  An  approximate  summation  pro- 
cedure was  employed  to  calculate  the  stress  distribution.   Deflections  were 
calculated  on  an  assumption  of  a  linear  variation  of  strain  vdth  time.   De- 
flections were  also  calculated  using  a  creep  equation  of  the  tj^ipie  of  Sq,  7 
with  a  linear  function  of  stress  for  e   and  an  exponential  function  of  stress 
for  m.  No  experimental  verification  was  reported,  however. 

In  a  paper  by  Marin,  Pao,  and  Cuff  (4),  a  method  was  described  for 
calculating  bending  creep  rates  and  deflections  in  a  beam  of  material  for 
which  the  creep  behavior  is  different  in  tension  and  compression.   The 
results  were  compared  with  data  on  creep  of  Lucite  and  Plexiglas,  Again  a 
linear  time  function  was  assumed  which  was  not  a  close  approximation  to  the 
observed  behavior.   Different  power  functions  of  stress  were  employed  to 
describe  the  elastic  strains  and  creep  strains.   In  the  analysis  the  elastic 
and  creep  strains  were  treated  separately.   The  analysis • predicted  a  differ- 
ent stress  distribution  and  position  of  the  neutral  axis  for  the  elastic 
component  of  strains  than  for  the  creep  (time  dependent)  component  of 
strains  in  the  beam.   Since  these  t^vo  stress  distributions  m.ust  exist 
simultaneously  it  does  not  seem  possible  for  them  to  be  different. 

In  a  paper  xvhich  had  just  become  available  at  the  time  this  manuscript 
was  prepared,  Pao  and  Marin  (l4)  presented  another  analysis  for  bending 


4. 
creep  of  Plexiglas  in  which  a  complex  equation  for  creep  of  the  form  of 
Eq.  6  was  employed.   This  equation  did  not  describe  the  initial  part  of 
the  creep  curves  much  more  closely  than  a  linear  equation.   Identical  creep 
relations  in  tension  and  compression  were  assumed  although  the  authors 
previous  data  on  Plexiglas  showed  this  to  be  untrue.   The  analysis  permits 
the  stress  distribution  to  vary  with  time  but  neglects  the  resulting  changes 
in  elastic  strains.   Other  simplifications  were  fovmd  necessary  in  order  to 
solve  the  problem.   The  resulting  equation  predicted  a  constant  creep  rate 
after  about  1000  hours  \vhich  did  not  agree  v^ith  the  data^  which  were  pre- 
sented to  11,000  hours,  in  which  the  creep  rate  decreased  continuously. 
The  paper  also  described  m.ethods  for  determining  creep  deflections  of  beams 
in  which  the  bending  moment  varied  along  the  beam, 

DERIVATION  FOR  BENDING  MTH  EQUAL  CREEF  IN  TENSION  AHD  COMPRESSION 
AND  EQUAL  TIME  DSi  SNDENT  AND  INDEPENPl^JT  CONSTANTS 

In  previous  studies  (9,  11,  15)  the  creep  curves  obtained  from  a  canvas 
laminate  tested  at  different  constant  stresses  were  found  to  be  closely 
described  by  an  equation  of  the  form 

e  =  e  +  mt^  =  (e  t  +  m't^)  sinh  ~  Eq,  10 

o         ^  o  a  ^ 

o 

xvhere  the  symbols  are  the  same  as  described  above.   This  relation  has  been 

shovm  (9)  to  be  related  to  the  activation  energy  theory  of  creep.   It  has 

also  been  shown  previously  (15)  that  an  equation  for  the  stress-strain  curve 

obtained  from,  a  constant  strain  rate  tension  test  could  be  derived  from  Eq,  10 

and  that  the  resulting  equation  was  in  close  agreement  Tirith  the  experimental 

data. 
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Fortunately,  bending  creep  data  were  also  obtained  for  material  of 
identical  composition  made  by  the  same  manufacturer  at  about  the  same  time. 
These  data  were  obtained  by  Marin  (20 ).  YiJhile  Marin  also  obtained  tension 
creep  curves  on  this  material  xvhich  agreed  in  general  with  the  tension 
creep  curves  in  reference  (9),  they  did  not  have  sufficient  self  consistency 
to  permit  evaluating  the  constants  in  Eq.  10. 

Dr.  Marin  in  a  private  communcation  provided  compression  creep  curves 
of  a  canvas  laminate  which  shov;  a  greater  creep  for  a  given  stress  than  the 
tension  creep  data.   However,  the  instantaneous  strains  were  also  50  percent 
great er^  ifnereas  the  stress-strain  curves  of  this  material-  in  tension  and 
compression  were  nearly  the  same  for  the  stresses  involved  in  the  creep 
tests  (12),  In  view  of  this  obserx'-ation,  the  following  derivation  has  been 
prepared  on  the  assumption  that  creep  in  tension  and  compression  was  iden- 
tical for  the  canvas  laminate, 

Eq,  10  accurately  represents  the  creep  data  under  conditions  of  constant 
stress.   It  probably  does  not,  however,  accurately  describe  creep  under 
varying  stress.   Thus,  if  the  stress  distribution  in  the  beam  changed  during 
creep,  as  described  by  Popov  (l9),  some  inaccuracy  might  be  introduced  by 
the  use  of  this  equation  unless  the  material  obeys  the  time-hardening 
law  (21),  which  is  probably  not  the  case. 

To  determine  the  bending  creep  relation  the  conditions  of  equilibrium 
and  geom-etry  must  be  satisfied.   From  the  geometry  of  a  beam  in  pure  bending, 
or  one  in  TAhich  the  distortion  due  to  shear  is  negligible,  the  following 
well-knovm  geom.etrical  relation  applies  : 

e  =  -^  Eq,  11 

P 

where  e  is  the  strain  at  a  distance  y  measured  perpendicularly  from  the 
neutral  axis  and  p  is  the  radius  of  curvature  of  the  centroidal  axis.   This 


6. 
relation  is  equivalent  to  assuming  that  plane  sections  remain  plane  -  a 
fact  which  has  been  demonstrated  for  plastically  bent  beams  by  Hadai  (22) 
and  for  creep  of  lead  beams  by  MacCullough  (16),  Also,  consideration  of 
compatability  of  strains  in  the  different  "fibers"  of  a  beam  of  great  length, 
bent  by  a  constant  bending  moment  into  the  arc  of  a  circle,  suggests  that 
plane  sections  must  remain  plane  at  all  points  along  the  beam  which  are 
remote  from  the  ends. 

The  conditions  of  equilibrium  for  a  beam  are  expressed  by  the  following 
equations  when  the  accelerations  involved  in  the  creep  are  so  small  as  to 
be  negligible  -  which  is  usually  the  case: 


oy(dA)  = 


Eq.  12 


a(dA)  =  0 


Eq.  13 


where  A  is  the  cross-sectional  area  of  the  beam,  h  is  the  depth  of  the  beam, 
and  c  is  the  distance  from  the  neutral  axis  to  the  tension  face  of  the  beam. 

Substituting  Eq.  11  for  e  in  Eq,  10  and  solving  for  the  stress  a  the 
following  is  obtained: 


o  =  c     sinh 
o 


n^ 


p(e  »  +  m»t  ) 
c 


=  a  sinh  (Ny) 


where  the  symbol 


N  = 


n' 


p(e  '  +  m't  ) 
o 


Eq.  14 


Eq.  15 


has  been  introduced. 


7. 

For  a  beam  of   rectangular  cross-section  of  width  b,    (dA)   =  b(dy),    in 
Eq.    12  and  13.      Making  this   substitution  and  introducing  Eq.   14  for  a  in 
Eq.    15,    the  follo7/ing  results: 


(dy)   ^  he  j  sinh~"^(Wy)(dy)  Eq.    16 

Upon  integration,  assuming  identical  creep  relations  for  tension  and  co-ii- 

h 
pression,  Eq.  16  shows  that  c  =  -^^  i.e.,  the  neutral  aris  is  at  the  center 

of  the  beam  and  does  not  change  position  Vifith  time  under  load. 

Making  the  same  substitutions  in  Eq,  12  and  changing  the  limits  in 

accordance  with  the  previous  paragraph  yields: 


/c  /c 

M  =  b  /    oy(dy)  =  ba   [  y  sinir"^(Ny)  (dy)         Eq.  17 
J-c  J-o 

After  integrating  and  rearranging  terms  this  becomes 


2M    f^.   N-2 
=  (lic) 


c  bc^ 


2(Nc)^  +  1|   sinh~''"(Nc)  -  (Nc  ).  j  l+dlc  )^  f    Eq.  18 
J  ^        ) 


which  for  convenience  in  computations  may  be  written 

2M    2e  Cosh  20  -  Sinn  28 

— — .  as  ,  , 

o  bc^       Cosh  26  -  1 


where  Q   =  Sinh"  Nc • 

Eq.  18  shows  that  the  bending  moment  M  is  a  function  of  K  for  a  beam 

of  given  size  and  material.  Thus,  for  a  constant  bending  moment  M,  N  must 

be  a  constant.  Hence,  from  Eq.  15  it  is  observed  that  the  curvature  —  of 

P 

the  beam  under  a  constant  bending  moment  varies  with  time  in  accordance 
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with  the  same  time  function  as  strain  varies  in  a  tension  creep  test.   That  is. 


-  =  N(e  '  +  m't^)  Eq.  19 


P 

In  this  equation  W  is  the  function   of  bending  moment,    M,    given  by  Eq.    18. 

The  most  convenient  way  of  determining  K  from  Eq,    18  is  by  means   of  a 

2M 
diagram  of  versus  Nc,   v/hich  is    easily  constructed  from  the   second 

o  bc^ 
o 

form  of  Eq.  18  as  shoi/vn  in  Fig.  1, 

From  Eq.  14  the  stress  distribution 

1      c 
y  =  —  sinh  —  Eq.  20 

o 

is   obtained.      Since  W  is  a  function  of  the  bending  moment,    size  and  material 
of  the  beam  and  independent   of  time   (according  to  Eq,    18),    it  follov;s   that 
the  stress  distribution  does  not  change  tvith  time  during  creep  under  a 
constant  bending  moment. 

APPLIGATIOIT  TO  CREEP  OF  CAiWAS  LAMINATE  IN  HIRE  BENDING 

Deflections 

Numerical  values  of  the  constants  in  Eq,  10  which  were  determined  from 
tension  creep  data  (9)  are  as  follows: 

e  »  =  m'  =  0.001875  in  per  in. 
o  ^ 

n  =  0.1183 
a     =  4000  psi. 

The  values  differ  from  the  values  determined  for  this  laminate  in  a  previous 
study  (15)  in  which  m'  was  0,0029  and  a  was  5230»   The  agreement  betT;een 
the  tension  creep  test  data  and  Eq.  10  using  the  present  constants  is  shown 
in  Fig.  2 


9. 

The  agreement  is  excellent  except  for  the  highest  stresses  of  5700  and 
6650  psi.   Since  the  specimen  at  the  latter  stress  fractured  at  550  hours, 
it  is  possible  that  the  fracturing  mechanism  may  have  influenced  the  creep 
results  at  the  highest  stress.   It  is  also  always  possible  that  specimens 
may  have  been  slightly  overstrained  during  loading,  although  great  care  vifas 
exercised  in  this  respect.   It  has  been  shown  (23)  that  such  overstraining 
alters  the  shape  of  the  creep  curve.   The  short  horizontal  lines  at  the 
left  border  of  Fig.  Z   indicate  the  first  test  points  for  each  test. 

To  calculate  the  deflection-time  curves  for  a  beam  of  the  dimensions 
employed  by  Marin  (20),  values  of  H  were  determined  from  Fig.  1,  correspond- 
ing to  the  desired  bending  moments.   These  values  were  substituted  into  Eq. 

19  together  virith  the  values  of  the  constants  given  above  and  the  curvature 

1 

--  was  determined  for  various  values   of  time  t.      Then  the  deflection  co  at 
P 

midspan  of  a  beam  of  length   i  subjected  to  a  uniform  bending  mo:?.ent  was 
calculated  for   each  curvature  from  the  well  Icnoira  geometrical   relation 

,/2 

6)  =  -^  Eq.    21 

8p 

The  resulting  theoretical  bending  creep  curves  are  shown  in  Fig.  3 
together  with  the  experimental  data  reported  by  Marin  (20)  for  a  grade  C 
canvas  laminate.   In  Fig.  5  the  solid  lines  show  the  bending  creep  pre- 
dicted by  Eq.  19  using  the  constants  employed  in  computing  the  solid  lines 
shown  in  Fig.  2,  The  black  dots  in  Fig.  3  represent  the  creep  that  would 
have  been  predicted  hj   using  the  constants  previously  em.ployed  to  fit  Eq. 
10  to  the  data  of  Fig.  2. 

The  short  horizontal  lines  at  the  left  border  of  Fig.  3  indicate  the 
deflection  predicted  at  zero  time  if  the  load  could  have  been  applied 
instantaneously  without  shock. 
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The  agreement  between  the  theory  given  by  the  solid  lines  of  Fig,  3 
and  experimental  data  is  reasonably  good.  Sone  of  the  discrepancies  which 
do  exist  may  be  explained  as  follows:  For  the  curve  at  a  bending  moment 
of  624  lb, -in,  the  experimental  data  indicate  a  constant  deflection  be- 
tween 200  and  900  hours.   Since  curves  for  higher  and  lower  momsnts  showed 
a  continuously  increasing  deflection,  it  seems  likely  that  experimental 
difficulties,  such  as  the  sticking  of  a  dial,  may  have  been  responsible  for 
the  constant  deflection. 

The  divergence  in  the  early  part  of  the  test  at  a  bending  moment  of 
767  lb, -in,  may  have  been  caused  by  overstraining  during  load  application. 
This  possible  explanation  seems  to  be  supported  by  the  fact  that  the 
difference  between  experimental  and  theoretical  data  at  this  load  decreases 
as  the  time  increases.   This  trend  is  in  accord,  with  the  observations  of 
Everett  (23)  and  MacCullough  (16)  on  the  effect  of  overstraining  while 
loading. 

The  experimental  data  for  bending  moments  of  917  lb, -in.  and  932  lb, -in. 
show  a  difference  that  seems  too  large  for  the  small  change  in  bending 
moment.   There  is  no  way  of  knowing  which  of  the  two  curves  is  m.ore  nearly 
correct. 


§tr 


esses 


The  stress  distribution  in  the  beam  under  a  bending  moment  of  679  lb, -in, 
v/as  computed  from  Eq.  20  and  is  shovai  in  Fig.  4,   The  maximum  stresses  in 
the  beams  were  also  computed  for  the  several  beams  for  which  data  are 
reported  in  Fig.  3,   These  stresses,  in  order  of  increasing  bending  moments, 
were  4920,  5640,  6060,  6740,  7940,  8,000  psi,  respectively.   It  was  observed 
that  the  two  highest  stresses  are  very  near  the  stress  in  the  tension  test, 
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Fig.    5,    at  T;hich  the   strain  increased  rapidly  with  small  increases   in  stress. 
This  may  account  in  par':  for  the  large  difference  between  the  creep  of  the 
tViTo  beams  having  the  highest  loads,    see  Fig.    3. 

RELATION   TO  TENSION  TEST  OF  CAIIVAS  LM^illvTATE 

In  a  prior  paper  (15)  an  equation  for  the  stress-strain  relation  in  a 
tension  test  v/as  derived  from  a  creep  equation  of  the  form  of  Eq.  10.  The 
resulting  stress-strain  equation  vms: 


e  =  e     smh  — 
o  o 

o 


1  + 


/  e  n 
/     o 


sinh 


et 


a 


1-n 


Eq.   22 


where  e  is  the  strain  rate. 

Excellent  agreement  was  found  betv\reen  the  stress-strain  curve  computed 
from  the  creep  data  for  the  canvas  laminate  used  in  the  above  study  and  the 
data  from  a  tension  test.   Since  the  constants  in  the  creep  equation  as 
employed  in  the  present  analysis  differ  from  those  used  in  the  previous 
paper,  the  stress-strain  relation  has  been  recomputed  from  Eq.  22  using 
the  new  constants.   The  theoretical  and  experimental  curves  are  compared  in 
F'ig,  5.   The  agreement  is  close  but  not  as  satisfactory  as  for  the  previous 
set  of  constants  (l5), 

noh-ibuform  bending  moment 

If  the  bending  moment  M  in  the  beam  varies  v/ith  distance  x  along  the 
beam,  then  M  may  be  expressed  in  terms  of  x  as  M  =  f (x).  Hence  from  Eq,  18 
N  is  a  function  of  x,  N  =  g(x).   If  in  Eq.  19  g(x)  is  substituted  for  N  and 


dfco 
dx^ 


,d"-w 


I  *  (i)^' 


-r  dx^  (for  small  deflections; 


;q,   23 


ISNSyEBSlW  OF  ILLW03S 


12, 
is  substituted  for  the  curvature  the  differential  equation 

l5e  =  g(x)  (e  •  +  m't'')  Eq.  24 

dx®         ° 

is  obtained  which  majr  be  solved  for  the  deflection  of  the  beam. 

Example  1,  Let  g(x)  =  Qx  where  Q  is  a  constant.   This  corresponds  to 

a  beam  having  a  bending  moment  M  v;hich  varies  with  distance  ::  according  to 

the  equation  obtained  from  Eq,  IS  by  substituting  Qx  for  K.   The  shape  of 

this  bending  moment  diagram  is  the  same  as  the  curve  shovm  in  Fig.  1  since 

the  ordinate  is  proportional  to  the  bending  moment  and  the  abscissa  is 

proportional  to  x.  Within  the  range  of  Nc  shown.  Fig.  1  has  approximately 

the  shape  bending  moment  diagram  v;hich  would  be  produced  in  a  cantilever 

beam  loaded  with  its  ovvn  weight  and  having  an  upward  force  at  the  free  end. 

At  a  given  time  t-^  the  deflection  of  the  beam  may  be  determined  in  a 

manner  identical  to  that  employed  for  elastic  beams.  From  Eq,  24 


—   =  Qx(e  •  +  m't,^)  =  QT  X  Eq.  £5 

dx^      ° 


where  T  =  s  '  +  m'tn^ 
Integrating  once 


^=QT|!.0 


If  the  beam  is  rigidly  clamped  at  x  =  L,  then 


dO) 

-—  -  0 

dx 


at  X  =  L 

and  C  =  -  QT  ~- 


Integrating  again 


x^      L^ 
w  =  QT  —  -  QT  J-  X  +  Ci 


15. 

But  <o  =  0  at  X  =  L  so  Ci   =  ^^  ,   hence,    co  =  QT  ~  -  QTL^^  |  +  ^~- 

which  is  the  deflection  curve  of  the  beam  for  small  deflections.     At  x  =  0 

the  deflection  co  is 

CO  =  ^  =  S^:.  (e    .    +  ^jt^^)  Eq,   26 

3  3  o  ■^ 

When  the  bending  moment  varies  with  distance  x  then  the  stress 
distribution  will  vary  also  with  x  in  accordance  with  the  expression  ob- 
tained by  substituting  g(x)  for  H  in  Eq.  20. 

y  =  -rir  sinh  ^  Eq.  27 

Example  2,  Eq.  18  usually  cannot  be  solved  for  N  as  an  explicit 
function  of  x.   Thus,  Eq,  24  cannot  usually  be  evaluated  by  integration, 
A  somev.'hat  different  approach  is  described  belox^: 

Consider  a  cantilever  beam  of  length  L  subiected  to  a  concentrated  load 
P  at  the  end.   Let  the  depth  be  2c,  and  the  beam  width  b.   Measuring  distance 
X  along  the  beam  from  the  free  end,  the  bending  moment  M  is 

M  =  Px  Eq.  28 

Differentiating  with  respect  to  x 

dM  =  Pdx  Eq,  29 

Eq.    14  may  be  written  for  the  stress   a-j_  in  the  extreme  "fiber"   as 

Nc   =  sinh  -^  =   sinh  S  Eq.    30 

a 
o 


where  S  =  — ^  ,  and  S  is  a  function  of  x  since  K  is  a  function  of  x.   Ec.  IS 
c 
o 

can  be  rewritten  as  follows  by  substituting  Eq.  50  and  cancelling  terms 

o  bc^  r  -| 

M  =  -2 I  s(2  +  csch^  S)  -  coth  S|  Eq.  31 

L  J 
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Substituting  for  M  from  Eq.  28  and  differentiating,  the  following  is 
obtained 

.2 


d:-  = 


o  bo' 
o 


coth^  S  -  S  coth  S  csch^  S 


_ 


dS 


Eq.    32 


From  Jiq.    10  and  11,    the  following  majr  be  written  for  the   extreme  fiber 


,  (e    '   +  m't'^)   sinh  S 

p       c  c 


Eq.    33 


After  eliminating  —  bet-^een  Eq,  23  and  33,  the  following  integral  can  be 
formed  in  terms  of  S  and  x 


dec 
dx 


'(e  '  +  m't°)  sinh  S 


dx 


Eq.  34 


Substituting  dx  from  Eq.  32  into  Eq.  34  and  integrating,  an  equation  for 

the  slope,  -r—  ,  of  the  beam  is  obtained 
'^    dx 


dco 


a  be 


^  =  -4—  (e  '  +  m't^)(cosh  S  +  S  csch  S  +  C,) 
dx    P     o  ^ 


Eq.  35 


where  C^  is  a  constant  to  be  determined  from  boundary  conditions. 

The  deflection,  w,  of  the  beam  may  be  found  by  multiplying  Eq,  35  by 
dx,  substituting  for  dx  from  Eq,  22  and  forming  the  integral.  After 
expanding  and  collecting  termn,  the  expression  becomes 

o  %%^  f  , 

CO  =  ~ (e    «   +  m»t^)    /  (coth^  S  cosh  S  -  S^  coth  S  csch     S 

p2      o  j 


+  Ci  coth^  S  -  Ci  G  coth  S  csch^  S)  dS 


Eq.  56 


The  indicated  integration  may  be  evaluated  with  the  following  result 
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c  %%'^ 


w  =-2 ^  (e   •   +  mH'') 


p2 


•   V,  Q  ^           S^         4.     S  cosh  S 
sinh  S  +  « = — -  + 

3  sinh     S       3  sinh^  S 


2          ^  S  ^  1   "^^         (-1)^  2    (g^P-^  >1)  g  _2p-l 
3  sinh  S       3       3   /_ (2p  -  1}   I  p 

p  — 1 J  C  m  »  % 


^  2   smh  S        2   ^.^j^^a  3^ 


Eq.    37 


where  S  =  —3'  ,    On    is  the  stress   on  the  extreme  fiber  of  the  beam  at   the 

o 

position  for  which  the  deflection  co  is  to  be  determined,   B     are  Bernoulli 
^  p 

numbers,   and  C^^  and  Cg  are  constants  to  be  determined  from  boundary  con- 
ditions.     The  series  appearing  in  Eq,    37  converges  for  S^<jr    • 

For  the  cantilever  beam  of  length  L  the  slope  is   zero  at  x  *  L,      Thus 
from  Eq.    35 

Ct    =  -  cosh  S^    -  S  csch  S^  Eq.    38 

Li  Li 

where  S     is  the  value  of  S  at  the  position  x  =  L.      The  value  of   S     may 
be  determined  from  Eq,    31  by  substituting  PL  for  M  and  solving  for  S     by 
trial#in  a  numerical  example. 

Since  the  deflection  of  the  cantilever  beam  is   zero  at  x  =  L,   the 
value     of  Cg  may  be  determined  from  t-q.    37  by  setting  co  =  0  and  S   equal 
the  same  value  S     as   above, 

Li 

To  obtain  the  deflection  co     at  location  x  along  the  beam,    substitute 
the  bending  moment  Px  for  M  in  Eq,    31  and  solve  for  the  corresponding  value 
of  S,    S   ,   by  trial.      Substitute  S      in  Eq,    37  together  with  C^  and  Cg  and 
compute  the  deflection  co  , 

It  can  be  shewn  by  application  of  L'Hospital's   rule  to  Eq,    31  that 

S  =  0  when  M  =  Px  =  0.     Thus  the  deflection,   o)  ,    of  the  free   end  of  the 

o 

beam  may  be  obtained  by  substituting  S  =  0  in  Eq.  37 
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%^-'^ 


O      DC 

(0     =  -2." (e    '    +  m't")  Co  Eq.    39 

o  p2        '    o  '      ^ 

The  above  solutions  are  valid  only  v;hen  S    <7r      as  noted  above. 

For  values   of     S     greater  than  about  5  the  hyperbolic  functions  may 

be  approximated  by  exponential  functions   to  obtain  a   solution.      Eq.    31   then 

becomes  approximately 

a  bc^ 
M  =  -~—  (2S  -  1)  So.   40 

Differentiating  and  substituting  Eq.  29 

dM  =  Pdx  =  a  bc^dS  So.  41 

o 

Svibstituting  Eq.    41  into  Eq.    34  and  integrating^   the   equation  for  slope 
of   the  beam  is   obtained 

~  =  ^4r  (e^'    -^  m't^)(e^   f  C,)  Eq.    42 

The  deflection  equation  is  obtained  by  substituting  for  dx  from  Eq.  41 

into  Eq,  42  and  integrating 

s-u  2„  3 
a  0  c 

0)  =  -2 (e  »  +  m»t^)(e  +  G^S  +  Co)  Ea.  43 

2P^ 

where  the  s3,Tabols  are  the  same  as  before. 

This  solution  also  involves  determining  S  as  a  function  of  x  through 

Eq.  28  and  40.   The  above  expression  will  be  a  good  approximation  to  Eq.  27 

when  the  value  of  S  is  greater  than  5  throughout  most  of  the  depth  of  the 

beam  and  throughout  most  of  the  length  of  the  beam.   Inaccuracies  will 

always  be  found  in  the  vicinity  of  the  neutral  axis  and  in  a  region  near 

zero  bending  moment  as  at  the  end  of  the  cantilever  beam. 
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EFFECT  OF  LIMITII-IG  VALUES  OF  THB  CREEP  CONSTAETS 

Examination  of  Eq.  10  shoivs  that  as  n  approaches  unity  the  creep  curve 

approaches  nearer  to  a  linear  curve  (constant  creep  rate)  v-hich  is  the 

behavior  expected  from  a  pureljr  viscous  material.  Also  as  a  becomes 

laree  so  that  —  is  small,  the  stress  function  becomes  nearer  a  linear 

o 
relation.   This  again  is  v;hat  would  be  expected  from  viscous  behavior. 

As  n  becomes  smaller  the  resulting  creep  curve  tends  to  rise  sharply 

and  then  level  off  to  an  apparently  constant  value  -  ^vhich  is  not  actually 

constant.   For  small  values  of  o  ,  such  that  —  is  larce,  the  hvperbolic 

o'  a        ' 

o 

sine  term  may  be  replaced  by   an  exponential  one. 

It  should  be  pointed  out  that  erroneous  results  are  obtained  by  sub- 
stituting — -  for  sinh  —  in  the  expressions  calculated  for  bending  creep, 

o  o 

Eq,  14  and  18,   This  does  not  indicate  an  error  in  the  equations  but  simply 

indicates  that  the  simplif j/-ing  assumption  for  small  stresses  cannot  be 
utilized  in  equations  v/hich  include  a  difference  of  t?ro  terms  of  like 
magnitude,  as  the  error  in  making  the  simplification  may  be  larger  than  the 
value  of  the  difference  of  the  t^vo  terms. 


18. 

REFERENCES 

1,  R,  W,  Bailey,  "TTtilization  of  Creep  Test  Data  in  Engineering  Design", 
Inst,  of  Mechanical  Engineers,  ¥ol.  131,  1935,  p,  131, 

2,  H.  J.  Tapsell  and  A.  E.  Johnson,  "An  Investigation  of  the  Nat'are  of 
Creep  Under  '^tresses  Produced  by  Pure  Flexure",  Inst,  of  Metals  Journal, 
Vol.  57,  Aug.  1935. 

3,  C.  C.  Davenport,  "Correlation  of  Creep  and  Relaxation  Properties  of 
Copper",  ASME  Trans.,  Vol.  60,  1838,  p.  A-55, 

4,  J.  iviarin,  Y„  H.  Pao,  and  G,  i^.  Cuff,  "Creep  Properties  of  Lucite  and 
Plexiglas  for  Tension,  Compression,  Bending,  and  Torsion",  Trans,  A31VIE, 
Vol.  73,  1951,  p.  705-719, 

5,  W,  Kauzmarja,  "Flow  of  Solid  Metals  from  the  Standpoint  of  the  Chemical 
Rate  Theory",  Trans.  American  Inst,  of  Mining  and  Metallurgical 
Engineering,  Institute  of  Metals  Division,  Vol.  143,  1941,  p.  57. 

6,  R.  G,  Sturra,  C,  Dumont,  sjid  F,  i,/[,  Koivell,  "A  Method  of  Analyzing  Creep 
Data",  Journal  of  Applied  Mechanics,  Vol.  3,  Wo* -2,  June  1936,  p.  A62, 

7,  H.  Leaderman,  "Creep,  Elastic  Hysteresis,  and  Damping  in  Bakelite  under 
Torsion",  Journal  of  Applied  :.xechanios.  Vol.  61,  June  1939,  p.  A79, 

8,  D.  Telfair,  T,  S.  Carswell,  and  H.  K.  Wason,  "Creep  Properties  of 
Molded  Phenolic  Plastics",  Modern  Plastics,  Vol,  21,  p,  137,  Feb.,  1944 

9,  W.  H.  Findley,  "Creep  Characteristics  of  Plastics",  1944  SjTnposium  on 
Plastics,  ASTM,  1944,  p.  18. 

10,  E.  P.  Popov,  "Correlation  of  Tension  Creep  Tests  with  Relaxation  Tests", 
Journal  of  Applied  Mechanics,  Vol.  69,  June  1947,  p.  A135, 

11,  W.  W.  Findley,  C.  H.  Adams,  and  'L    J.  Worley,  "The  Effect  of  the  Creep 
of  Two  Laminated  Plastics  as  Interpreted  by  the  Hyperbolic  -  Sine  Law 
and  Activation  Energy  Theory",  Proceedings,  ASTM,  Vol.  48,  1948,  p.  1217. 

12,  W.  IJ.  Findley  and  W.  J.  Viforley,  "Mechanical  Properties  of  Five  Laminated 
Plastics",  MCA,  Tech.  Note  1560,  August  1948. 

13,  W.  J.  Worlej'-  and  W,  N,  Findley,  "The  Effect  of  Temperature  on  the  Creep 
and  i:^ecovery  of  a  Uelajrdne- Glass  Fabric  Laminate",  Proceedings,  ASTIJ, 
Vol,  50,  1950. 

14,  Y,  Pao  and  J.  Marin,  "Deflections  and  Stresses  in  Beams  Subjected  to 
Bending  and  Creep",  Applied  Mech.  Division,  ASIioE  Paper  l-Io.  5^-APL^34, 
Presented  June  1952, 

15,  W.  N.  Findley,  "Derivation  of  a  Stress-Strain  Equation  from  Creep  Data 
for  Plastics",  Proc.  First  U,  S.  riational  Congress  of  Applied  Mechanics, 
June  11-16,  1951,  p.  595. 


19. 


16.  G,   H,   MacCullough,    "An  Experimental  and  Analytical   Investigation  of 
Creep  in  Bending",    Trans.   ASllffi,    Vol.    55,    1933,   p.   APM  55-9-55. 

17.  J.   Marin  and  L,    E.    Zwissler,    "Creep  of  Aluminum  Subjected  to  Bending 
at  Normal  Temperatures",    Proc.  ASTII,   Vol.    40,    1940,    p.    S37. 

18.  J.   Iviarin  and  G.   Cuff,    "Creep-Time  Relations  for  Polystyrene  under 
Tension,   Bending,    and  Torsion",    Proceedings  ASTM,   Vol.   4S,    1949, 
p.    1158-1174. 

19.  E.   P.   Popov,    "Bending  of  Beams  v/ith  Creep",   Jour,   of  App.   Physics, 
\^ol.    20,    1949,    p.    251-256. 

20.  J.   Marin,    "Static  and  Dynamic   Creep  Properties   of  Laminated  Plastics 
for  Various  Types   of   Stress",   NAGA,    Tech.   Note  1105,   Feb.    1947. 

21.  I.   Roberts, "Prediction   of  Relaxation   of  I.Ietals  from  Creep  Data", 
Proceedings  ASTM,   Vol.   51,    1951^   p.   811o 

22.  A.   Nadai,    "Plasticity",   The  McGraw  Hill  Book  Co.,    Inc.,    1935,    p.    121. 


23.     F.   L.   Everett,    "Strength   of  Ivlaterials   Subjected  to  Shear  at  High 
Temperatures",    Trans.   AS'VIE,    Vol.    53,    1931,   p.    126. 


t!A 


(/) 


CO 

(O 

llJ 

q: 

u.- 

o 

h- 

C/l 

N- 

h- 

1- 

2 

< 

LlI 

tr 

Q 

UJ 

UJ 

Li_ 

H 

L- 

CO 

o 

UJ 

1- 

h- 

UJ 

< 

o      ^ 

2 

-     o 

S 

-      CA) 

< 

(U         2 

1 

E       UJ 

1-       1- 

CO 

< 

2 

> 

< 

cn 

o 

K 

co 

o 

UJ 

1 

H 

UJ 

Q 

CL 

< 

UJ 

q: 

UJ 

CD 

(r 

o 

or 

o 

(\J 

u. 

CI) 

IJL 

^        5!:!        o        00        to 
u!  jad  uigOl  'uiDJjs 


o 

si-" 

ID 

ro 

O 

in 

CM 

q 
cm' 

IMS 

ID 

d 


<NJ 


CO 

> 


00 


< 

O 
UJ 


S5     o 

u_      cc 


^ 

1^ 

0^ 

^ 

^ 

o 

^"^ 

^  -^ 

^T^ 

o 

O  ~ 

0) 

o 

•  1                                0 

o 

e    <!• 

o* 

o 

ll 

ll 

o 

() 

o 

'^ 

"^ 

o 

o 

. 

i 

ll 

QS^- 

CX) 

h- 

C/) 

< 

o 

o 

!? 

2 

> 

o^ 

o 

<> 

o 

H 

UJ 

z 
< 

o 

s 

b° 

b'' 

1^ 

UJ 
U. 

o 

C 

*^   K^        i^ 

?^  {^ 

j:^ 

<l> 

^ 

k 

o 

o^ 

D  r-^°     ^c 

«^     CVi 

of^ 

/^ 

^ 

s^ 

1 

C 

^    po 

p 
o 

^ 

Si 

0) 

5 

1  - 

CD 

Q 
< 

UJ 

o 

<  u: 

o 

(> 

o 

q:   o 

o 

v^ 

C3 

g- 

o 

e>   1^ 

o 

o 

<> 

o 

o 

1 

1- 

§. 

O 

z 

D 

^ 

<o 

<u 

o 

O    H 

n 

1        "^ 

/-V 

1 

lO 

•* 

E 
i- 

UJ 

m 

NTS    F 
TED    A 

0 

ii       ° 

o 

> 

3 

• 
o 

0 

o 

1 

o 

to 

0) 

o  _ 

';^ 

C/) 

UJ     (/) 

c 

5          11         ° 

p 

o 

o 

2    UJ 

11 

1       ^- 

fc 

s 

h- 

O    H- 

11 

o 

cn 

2 

o      H 

3 

o 

D 

^ 

UJ 

UJ 

11       ° 

o 

o 

• 

m^^\ 

1- 

CD     1— 

o        11       ° 

ol 

=    \ 

ro 

Q- 
UJ 

2    < 

11     ° 

ol 

'  1    \ 

Q     ? 

o       U      o 

° 

0 

UJ 

2     2 

o        U     o 

1             < 

\\    ' 

o 

UJ     < 

GQ     J 

o          \\ 

o  1 

CJ 

O           U       o 

•o\      • 

U  .1 

W    o 

o             W 

I           1 

ro 

\\° 

o    \ 

r  \        1 

— 

(D 

o        \\   o 

\  \          T 

°  \ 

o             o  \ 

\  \      L 

U- 

o          \ 

v°      A 

•W\ 

o 

^       -0^ 

^^4^ 

O 

.,        ,,   _                            .    .,.,.        ,    ,          ,1 

-ITT" 

lO 

o                     i^                    O                    " 

0 

c 

> 

(XJ 

<\J 

" 

"■ 

— 

uj^OI    'uoijoaiiaa 


o 

o 

Theory    from 
eauation  22 

X 

N 

^ 

X 

^ 

^^^o 

^          0 

o 

2.  o  ^ 

~--^, 

a> 


00 


CO  lO  <^ 

jsd    01    '  ssajis 


ro 


CO 


00 


l£) 


CO 


<s> 


CO 


d   O 


CO 
LU 
> 

q: 

Z) 

o 


< 
q: 


OJ 

a> 

d 

H 
1 

o 

to 

'o 

in 

« 

LU 

03 

c 
o 

q: 

1- 

CO 

o 

I- 
< 

(- 

CO 

ID 
CD 


o 

I- 
< 

Ijj 


CO 

< 
O 

O 
I 

LiJ 
Q 
< 

a: 
cr> 

cr 
o 


R 

"^ 

x^ 

isd 
)09      00 

N 

\ 

00      40 
Stress , 

\ 

§ 

lo      ^     to      Oj     ;:-     J 
Q     d    Q     a     ^     ^ 
ui  'uuoag  jo  duij 

\     ^      c 

A  a    c 

3      Q       C 

1 

^   \ 

8 

k 

O 

1 

\ 

^ 

^ 

oc 

o 

li- 

ar 

LU 
UJ 

or 

LU 

o 

o 

CD 

-z. 

CD 

(T 

2 

3 
Q 

Q 

2 

^ 

LU 
OQ 

<. 

LlJ 

1    ■  0° 

OQ 

t^V 

2  <  2 

2  — 

T' 

^  P 

O 

1  ^ 

h- 

3 

CO 

c 

OJ 

or 

§ 

o 

1- 

2 

V) 

(O 

< 

Q 

o 

1 

Q. 

E 

CO 

o 

Ci 

CO 

LU 

LU 

Q 

ir 

< 

1- 

CC 

C/) 

CD 

CD 


679 


£^ 


^r^--^ 


no. 2 


CONFE 


15niit»eeir*^* 


l^ibrarr 


iiLiilii. 


PRELIMINARY  CONSIDERATION  OF  LINEAR  AND 

NONLINEAR  MECHANICAL  MODELS  FOR 

CREEP  OF  PLASTICS 


^ 


by 
E.  J.  Scott  and  W.  N.  Findley 


A  Research  Project  of  the 

DEPARTMENT  OF  THEORETICAL  AND  APPLIED  MECHANICS 

UNIVERSITY  OF  ILLINOIS 


Sponsored  by 

PICATINNY  ARSENAL,  ORDNANCE  CORPS 
DEPARTMENT  OF  THE  ARMY 

Interim  Report  No.  2 

on 

The  Relationship  Between  Time  Sensitive 

Mechanical  Properties  of  Plastics 

Contract  No.  DA-n-022-ORD-401,  Project  No.  TB4-721 


Urbana,  Illinois 
January,  1954 


Return  this  book  on  or  before  the 
latest  Dat^j;^ 


University  of  Illinois  Library 


^- 


^-M^i^^^"^' 


^  1^^' 


Interim  Report  No .  2 
on  a  research  project  entitled 
Study  of  Relationships  Between 
TJme  Sensitive  Mechanical 
Properties  of  Plastics 

Project  Supervisor,  W.  N.  Findle-; 


PRELIMINARY  CONSIDERATION  ON  LINEAR  AND  NONLINEAR 
MECHANICAL  MODELS  FOR  CREEP  OF  PLASTICS 


E.  J.  Scott 

Associate  Professor  of  Applied 
Mechanics,  Kansas  State  College 
(formerly  University  of  Illinois 

W.  N.  Findley 

Research  Associate  Professor  of 
Theoretical  and  Applied  Mechanics 
University  of  Illinois 


Department  of  Theoretical  and  Applied  Mechanics 
University  of  Illinois 


THE  LIBRARY  OF  THf" 
SEP  7  1954 
UNIVERSITY  OF  iU'hUiS 


^7? 


SUMriARY 

A  preliminary  ana.lysis  is  presented  for  three  four- 
element  models:   (a)  linear  Maxwell  and  Voigt  units  in 
series  with  two  masses,  (b)  two  nonlinear  Maxwell  units  in 
parallel,  coupled  by  a  mass,  and  (e)  nonlinear  Maxwell  and 
Voigt  units  in  series  vjith  two  masses. 

The  equations  of  m.otion  were  derived  in  each  case  and 
the  solutions  of  the  resulting  differential  equations  were 
indicated.   Complete  solutions  for  particular  problems  have 
not  been  included. 


CONCLUSIONS 

The  solution  of  mathematical  expressions  describing 
the  behavior  of  mechanical  models  composed  of  nonlinear 
elements  is  possible,  although  involved.   Since  these  models 
may  describe  the  behavior  of  plastics  more  closely  than 
models  composed  of  linear  elements  they  may  permit  more 
accurate  prediction  of  the  time  sensitive  behavior  of  plas- 
tics under  different  loading  conditioas.  Hence,  further 
study  of  such  models  seems  warranted  and  a  comparison  of 
the  predicted  behaviors  with  actual  material  performance  is 
desirable. 
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Various  types  of  mechanical  models  have  been  proposed  in 
the  past  for  describing  the  time-sensitive  behavior  of  material* 
These  have  be6n  composed  of  various  combinations  of  linear 
springs  and  linear  dash  pots,  that  is,  viscous  elements  whose 
resistance  is  proportional  to  the  velocity  of  separation  im- 
posed on  the  dash  pot.  Recently  nonlinear  dash  pots  were  studied. 

These  models  admittedly  do  not  in  any  way  represent  the 
actual  construction  of  molecular  materials.   In  fact,  models 
previously  considered  neglect  the  mass  of  the  material.  The 
mechanical  behavior  of  these  models,  however,  can  be  made  to 
approximate  some  of  the  observed  behavior  of  polymers  and  in 
fact  analysis  of  these  models  has  permitted  prediction  of 
mechanical  behavior  of  polymers .  The  use  of  these  models  con- 
sisting of  elastic  and  viscous  elements  has  led  to  the  descrip- 
tion of  a  class  of  material  having  characteristics  describable 
by  these  models  as  "visco  -  elastic"  materials. 

It  has  been  found  difficult  to  accurately  describe  behavior 
of  many  materials  by  a  single  element  composed  of  a  small  num- 
ber of  units .  So  a  large  number  of  elements  identical  in  form 
but  having  different  constants  has  often  been  postulated.  This 
has  the  effect  of  making  possible  a  distribution  of  "relaxation 
times". 

It  may  be  that  the  multiplicity  of  linear  elements  may 
be  replaceable  by  a  single  non-linear  element.  For  this 
reason  the  following  studies  have  been  undertaken. 

In  the  present  report  a  preliminary  analysis  is  presented 
for  three  four-element  models:   (a)  linear  Maxwell  and  Voigt 
units  in  series  with  two  masses,  (b)  two  non-linear  Maxwell 


-  2  - 

units  in  parallel,  coupled  by  a  mass,  and  (c)  non-linear 
Maxwell  and  Voigt  units  in  series  with  two  masses. 
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Linear  Maxvjell  and  Voigt  Units  in  Series  v^lth  Two  Masses 

A  model  consisting  of  a  Maxwell  unit  (spring  and  dash- 
pot  in  series)  and  a  Voigt  unit  (spring  and  dashpot  in  parallel) 
in  series  with  a  mass  m   at  either  end  of  the  Voigt  unit  is 
diagrammed  in  Pig.  1.  This  type  model  is  particularly  inter- 
esting since  it  predicts  in  a  general  way  not  only  creep  but 
recovery  on  removal  of  load. 

In  this  section  the  differential  equations  of  motion  of 
this  model  are  derived  and  solved  for  the  case  of  linear  springs, 
k,  and  dashpots,  c.   Tv/o  forcing  functions  F(t)  of  practical 
interest  are  considered. 

From  the  free  "body  diagrams  of  three  portions  of  the  model. 
the  following  equations  of  motion  may  be  written: 

\^^X^{t)   =  o^^i^it)    -  X^(t)]  ,  (l) 

^X^{t)   =  kjx^(t)  -  X2(t)j  +  C2[x^(t)  -  i^it)], 

-  c^[x2(t)  -  X^{t)^   ,  (2) 

X^(t)  =  F(t)  -  kJx^(t)  -  X2(t)]  -  C2[x^(t)  -  X^Ct)]  (3) 


m 


nig 


Solving  equation  (l)  for  X-,(t)in  terms  of  X2(t),  we 

obtain  _^^  ^   ^^ 

c   ^         c 
X^(t  =  e  ^  f   e  ■'■  X2(T)dT  .  (4) 


Also,  using  integration  by  parts,  we  have 

"^     -iT  -it  ,  r    — T 

c   .  c  k  I   c 

e  ^  X2(T)dT  =  X2(t)e  ^  -  X2(0)  -  J:   e  ^  X2(T)dT.   (5) 
o  Vo 
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Therefore,  from  (4)  and  (5) 
X^(t)  =  X^Ct)  -  e  1   X2(0)  -  J^  e  ^^   e^X2(T)dT.   (6) 
Now,  from  (l)  and  (6),  we  have 

X2(t)-  X^(t)  =  .^X^Ct)  =  ^X^Ct)  -  ^  e   ^  X2(0) 

-  -^  /   e   ^       X2(T)dT.  (7) 

o 

Substituting  this  result  in  (2),  we  get 

it 

•  .  c 

rn-j^Xg  =  kgCX^  -  Xg)  +  c^CX-j  -  Xg)  -  k-^Xg  +  k^^e  ^  X2(0) 


k'  r   -^^'  -  "^ 

+  e^/   e  ^       X2(T)dT.  (8) 


We  now  make  use  of  the  Laplace  transform  method.  To 
that  end,  let  L  rx^(t)]   =  x^  (p)  ,  L  [ XgCt)]   =  x^{p)   and 
LfF(t)]   =  f(p).  Then  the  Laplace  transforms  of  (8)  and  (3) 


are 

m 


if  P%(p)  -  pXgCo)  -  XgCo)   =  kgFx^Cp)  -  X2(p)j 
+  Cglp  x^(p)  -  X^(0)  -  pXgCp)  +  XgCo)!-  k^^XgCp) 

k3^X2(0)     k^  x^M 

"*'  ^ —  "*■  c~  E —  '  (^^ 

p  +-i        -^   p  +  ^ 
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and 


m. 


[^1 


[p%(p)  -  PX3(0)  -  X^(0)J  =  f(p)  -  k^fx^Cp)  -  X2(p)j 

1  (10) 

px^(p)  -  X^(0)  -  px2(p)  +  XgCO)  I  . 

Rearranging  the  terms  of  equations  (9)  and(lO),  we  have 

2 


p  +  CoP  +  (kn  +  ko)  - 


^1 


XgCp)  -  (CgP  +  k2)x^(p)  = 


k^^XgCO) 


ni^pX2(0)  +  m-j^X2(0)  -  0^X^(0)  +  02X3(0)  +  ^  ^^  ,  (ll) 


P  +  c- 
^1 


-(C2P  +  k2)x2(p)  +1  ni2P' 


+  C2P  +  kg 


]  X3(p)  = 


in2PX^(0) 


+  m2X^(0)  +  03X^(0)  -  02X2(0)  +  f(p).   (12) 


Equation  (12)  when  solved  for  x-,(p)  yields  the  result 

(C2P+  k2)x2(p)     m2PX^(0)+  m2X-5(0)+  C2X  (O) -02X3(0 )+f(p) 
(p) 5 +  5 


^  mgP  +  C2P  +  k2 


m^-p     +   CoP  +  K 


(13) 


This  result  when  substituted  in  (ll)  and  simplified  gives: 

(c^p  +  k^)(c2P  +  k2)(m2PX^(0)+  m2X^(0)  -  03X2(0)+  f(p)) 
^  pD(p) 

(c^^p  +  k^)(m2P^+  02P  +  k2)(m^pX2(0)+  m^X2(0)-^^^0^°2^'^) ) 

pD(p) 


C3^k^X2(0)(m2P  +  O2P  +  kg) 
pD(p) 


(1^) 


-6- 
where 


r      h  "5 

o 

+  (^2^1^!  "*"  "^2^1^2  '*'  "'2°1^2  ■*"  "'l^l^2  "^  '"l^A^P 

+  (k^k  m  +  k^c^Cg  +  m^k^k^)?  +  c^k^k^l  .  (15) 


Substituting  this  result  in  (13)^  we  have 
X  (p)  =  (c3^P+k^)(c2P+k2)^(m2PX^(0)+ni2X3(0)+C2X^(0)-C2X2(0)+f(p)) 

p(m2P  +  CgP  +  k2)D(p) 

(c^p+k^)(c2P+k2)(m-j^pX2(0)+m^X2(0)-C2X^(0)+C2X2(0)  ) 

pD(p) 


+  pI5Tpl + 5 ; 

^^^^  nigP  +  C2P  +  k2 

(16) 
Making  use  of  the  Complex  Inversion  Integral  we  find  that 
X2(t)  = 

1  ^^'^^'^(c^p+k^^ )(c2P+k2 )^pX^(0)+m2X^(0)+C2X^(0)-C2X2(0)+f (p))  ^^ 


27ri|  pBTpT 

c-i«?D 


e^  dp 


1  r°"*'^'^(  c-^p+k^^ )  (m2P^+C2P+k2 )  (m^pX2(0>f  1113^X2(0) -02X3(0) +02)^(0))  . 
-^c-ioo 
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.c+io«?  2 

1  r       (c^p+k^)(c2P+k2r(m2PX3(0)+m2X^(0)+C2X5(0).C2X2(O)+f(p))  p^ 

^ pdTigP"^  +  CgP  +  kgJMp) ^  ^P 

^  J_r   (c^p+k^)(c2P+k2)(m^pX2(0)+m^X2(0)-C2X5(0)+C2X2(0))  pt 
J  pD(p) 

1   /^^"^^^c^k^X2(0)(c2P+k2)  pt 

+  ''2¥r-j       pi^(F) ^  ^'^ 


i2PX^(0)+m2) 
c-i<^  rngP"-  +  CgP  +  kg 


+  if    ni2PX3(0)+m2X3(0)  +  C2X^(0).C2X2(0)+f  (p)  ^ 

Siriy    ^ ^ — 2 i e^  dp     (lo; 


The  form  of  the  solution  will  depend  on  the  function  f (p) 

2 
and  on  the  nature  of  the  roots  of  D(p)  =  0  and  nigp  +  C2P  +  k2  =  0 

Once  these  are  known  the  residues  of  the  above  integrands  may 

be  computed  at  those  roots  which  are  the  poles  of  the  integrands. 

The  nature  of  the  roots  of  D(p)  =  0  can  be  determined  by 

studying  its  discriminant  which  is 

A  =  I^  -  27J^,  (19) 

where 

T   .2n,  »  „  V   ('"l"'2''l+°l'^2'"2+°l°2"'l)('^lV2^'^l°l°2-"l''A) 
I  =  c^m-j^mgk^kg jr 

o 

(m2k^c^  +  m2k-j^C2  +m2C^k2  +  "^-^^^^   "^  "'l^^A^ 
+   1^ '         (20) 


■K4 
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and 

2 

+  ^gCm-j^nigk^+c-j^Cgmg+c^Cgni^ )  (k-^^k^m^+k^c^c^-pm-j^k^k^ )  (m2k^c-j^-HTi2k-j^C2 

1  2 

+  in2C-j^k2  +  '"i^i^2"*'"^1^1^2^"  j^c^mT^ni2(k^k2ni2+k-j_c^C2+m^k^k2) 

-  ■^g-(n]2k^c^+m2k^C2+m2C^k2+ni^k^C2+m^c^k2) 

-  Ar^  (m^m2k3_^-C3_C2m2+c^C2ni^)^  (2l) 

We  novj  have  the  following  cases: 

(1)  If  two  roots  are  equal  A  =  0, 

(2)  If  A  <  0,  two  roots  are  real  and  two  Imaginary 

(3)  If  A  >  0,  the  roots  are  all  real  or  all  Imaginary. 
In  order  to  apply  this  solution  to  creep,  F(t)  is  made 

equal  to  a  constant  F  for  t  >  0  and  X^(t.).  is  -detemiiiied.  --Values 
of  the  constants  m,,  m^,   e^  ,  e^,  k,  and  k^.  must  be  selected  to 
fit  the  characteristics  of  the  material  under  cjonsideration.  .. 
Also  the  following  initial  conditions  must  obtain: 

X^(0)  =  XgCO)  =  X^(0)  =  V^^l  and  X2(0)=  X^CO)  =  0   (22) 

For  stress  relaxation  X,(t)  is  made  equal  to  a  constant 
X,Q  for  t  >  0  and  F(t)  is  determined  for  the  .sajne  values  of  the 
constants  and  the  following  initial  conditioxis:  X,(0)  =  X2(0)  = 
X^Q,  F(0)  =.  X^^k^,  X^iO)   =0. 

The  solution  has  not  y$t  been  evaluated  for  a  specific 
problem. 
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Two  Non-linear  Maxwell  Units  in  Parallel,  Coupled  by  a  Mass 
Figure  2  illustrates  the  model  considered  in  this 
section.  Models  such  as  this  composed  of  linear  units  have 
been  used  to  describe  creep  behavior.  A  model  of  this  type  is 
also  capable  of  representing  recovery.  This  may  be  explained 
as  follows.   If  the  springs  and  dashpots  have  unequal  constants, 
then  the  deflection  of  the  two  springs  after  a  time  Interval 
under  load  vjill  be  unequal.  Thus,  when  the  load  is  released 
one  spring  will  have  a  residual  tens ion- -the  other  a  residual 
compression.  These  residual  internal  forces  will  operate  on 
the  dashpots  so  as  to  relieve  these  forces  and  in  so  doing  will 
produce  the  motion  of  recovery. 

As  far  as  is  known  the  non-linear  model  of  Fig.  2  has  not 
been  considered.   Since  the  behavior  of  plastics  appears  to  be 
non-linear  in  character  it  may  be  that  properly  chosen  non-linear 
elements  will  permit  a  more  accurate  description  of  the  observed 
behavior  and  hence  a  more  accurate  prediction  of  behavior  under 
other  circumstances. 

In  a  previous  paper   a  stress -strain  equation  was  derived 
which  was  in  good  agreement  with  experimental  data  from  a  canvas 
laminate.  The  time -independent  portion  of  this  equation  was  a 
hyperbolic  sine  function.  Hence,  it  seems  reasonable  to  express 
the  deflection  X  of  the  non-linear  springs  as  hyperbolic  sine 
functions  of  force  f : 

X^  =  b^  sinh  f^/a^  for  spring  1,  (25) 

Xg  =  b^  sinh  fg/ag  for  spring  2.  (24) 

where  a^,  a^,  b  ,  b^  are  constants  of  the  springs. 
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Also  creep  data  on  plastics  have  been  found  to  be  in 

accord  with  the  stress  dependence  predicted  by  the  activation 

2 
energy  theory  of  creep  .  This  stress  dependence  is  a  hyperbolic 

sine  of  stress.  Thus  it  seems  reasonable  to  express  the  velocity 

of  flow,  X  -  Xv)j  of  tl^s  time  sensitive  elements  of  the  model 

l^y  hyperbolic  sine  functions  of  force,  f : 

X  -  X^  =  K^  sinh  f^/^^  for  dashpot  1         (25) 

X  -  Xg  =  Kg  sinh  fi^/Sg  for  dashpot  2,        (26) 

where  K-,,  K^,  6-,,  Sp  are  the  constants  of  the  dashpots . 

By  forming  the  inverse  of  equations  (23-26)  so  as  to 
solve  for  the  forces,  and  applying  these  forces  to  free  body 
diagrams  of  portions  of  the  model  the  following  equations  of 
motion  of  the  system  are  found: 

6-j^  sinh'"-^(X  -  X^)/K-j^  -  a^  sinh"-*-  i:^/X>^--0  (2?) 

Bg  sinh"-^(X  -  i^)/K^   -  a^  sinh""^  ^2''^2  ^  °       ^^^^ 
F(t)  -  3^  sinh'-'-(X  -  X^^)^^  -  6^   sinh"-^(X-X2)/K2 

•  m 

=  mX  (29) 

These  equations  will  be  solved  for  the  special  case  that 

the  springs  are  linear,  m  =  0  and  F(t)  =  F^  (a  constant)  i.e., 

for  constant  load  creep.   If  the  springs  are  linear  Eq.  (23)  and 

(24)  become 

^1  =  ^l^/^l  =  ^l/^l  (50) 

^2  "  ^2^2 '^^2   "  ^2 ^^2  (5^^ 
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where  k-,  =  a, /b-.  ,  kg  =  apA)^,  Equations  (27-29)  may  then  be 

written 

Q^   sinh"-^(X  -  X^) /K^   -   k^X-j^  =  0  (52) 

62  sinh"^(X  -  X2)/K2  -  k^Xg  =0  (53) 

Fq  -  ^1^1  -  ^2^2  =  ^  (^^^ 

Inverting  Eq.  (32)  and  (33)  and  substituting  Eq.  (3^)  for 
Xg  Eq.  (32)  and  (33)  become 

X  -  X-j^  =  K^  sinh  k^X^/p.-j^  (35) 

X  -  Xg  +  Kg  sinh  (F^  -  k-j^X^)/s.  (36) 

Differentiating  Eq.  (3^)  yields 

^A  "^  ^2^2  "  ^  ^^'''^ 

Substituting  Eq.  (37)  in  Eq.  (36)  and  eliminating  X  between 
Eq.  (35)  and  (36)  gives 

K  k  K  k 

^1  =  Ep%  ^i"h  (P„-  k^Xi)/S,  -  jA^  sinh  k^X/B^      (38) 

Expanding  the  first  hyperbolic   sine   in  Eq.    (38)   and  collecting 
terms 

^1  =  4^  '■^^^  ^o'^\  °°^h  k^X^/6;L— -f-(Ki*K2CoshP/ei>lnhlft/9^ 

(39) 
or 

^1   ~     ^1  ^°^^  ^'I'^l  '*'  ^^2  ^^^^  ^1^1  ^^^^ 


LIBRARY 

uNivERsmr  OF  wm^ 


'■.-.. I      l 


.,  ;■  ■  s  I.  r-n 


M-V, 


'■4.7 


aiina-; 


-J  J 


.  X  ...'■ 


ir 
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and  a-,  =  k-,/6-,  . 

In  order  to  solve  the  non-linear  equation,  Eq.  (4o)  the 
reversion  method  may  be  employed.   Then  after  solving  Eq.  (4o), 
Xn  and  X^  as  functions  of  time  may  be  substituted  in  Eq.  (55) 
and  X  may  be  determined  by  integration. 

The  reversion  method  requires  the  series  expansion  of  the 
functions  on  the  right  hand  side  of  the  equality  sign  of  Eq.  (^O) 

as  follows:  ^  ^  u         55   . 

a  X^   °^1^1   °^i-^l 
11^   cosh  a-j^X^  =  jjL^(l  +  g/  +  -j^r-   +  -57-  +  .  .  .  ) 

^"^  a^X^   a5x5   aM 

lig  sinh  a^X^  =  \i^{a^Y.^   +  -jj—  +  -^j-  +  "yt-  +...). 

Consequently  substituting  in  Eq.  (^O)  yields 
^1  =  ^1  +  ^a^l^l  +  5T  »?i  +  7T  »l^i  +  TTT  "l^l  +  3T  "Pi 

+  ST  4^f  +  BT  "fxf  ^  ?T  "14  +  •  •  •  '  ("^l) 

or,  for  the  purpose  of  the  method  at  hand. 

The  reversion  method  obtains  closer  and  closer  approxima- 
tions to  the  solution  of  the  non-linear  equation  by  solving  the 
following  series  of  linear  equations 

^^  '   ^2^1)^11  =  ^^1,  (^^) 


mlv 


Xrtm:.-, 


-r.;.'.'   J- 


I  r  n.  \ 


\  :■  :• 


y-a    •■■:'B->-' 


;  r- 


.<<:t    '£• 


':.  rc 


'p:y 
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d  ^1  ^  ^   '^2  5  2       ^1  2  /  2  \ 

(^  -  M.2ai)X-L4  =  i|-r  aiXii+  2T^l%1^12  +  2T  "l^^^ll^l^"^  ^12^' 

etc.  (^6) 

where  the  solution  to  Eq.  (42)  is  given  by 

X^(t)  =  X^3^(t)  -  X-^2^^)  -^  X^^(t)  +  X-^^(t)  +  .  .  .     (47) 

The  work  involved  in  this  method  is  great  so  the  present 
example  will  be  carried  on  far  enough  so  that  the  procedure  is 
clear.  Further  approximations  can  be  found  by  using  additional 
formulas  given  in  the  reference' . 


The  solution  of  Eq.(4;))  is 

tigO'i^   1^1 

'-2'"1 


^ii(^)  ==  -1- ' '    ~  nk  ^  ^''^ 


If  the  initial  conditions  are  X-,(0)  =  X,-,(0)  =0.,  then 


^1 

c,  = and  therefore 

1    M-2°'l 

X3^^(t)  =  ^L^d  -  e^"4^),  (49) 

where 

11^  =  V^^a^     and  n^  ^  ^^;^^^ 

Equation  (49)  constitutes  the  first  approximation  to  the 

solution  of  Eq.  (4o).  The  second  approximation  is  obtained  by 

solving  the  equation 

2 

(-dT  -  f^4)^12^^)  =  -2T^  ^^1(1  -  e^^^^i^  .  (50) 
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which  is  obtained  from  equation  (44).  The  solution  of  Eq.  (50)  is 


X;^2(^)  = 


1x3  a^  2.1 


-  -2~r  '^^-^^-r  "^  ^^^ 


^4^  .  l^2^,t)  ^   /4^ 
M-4         2 


If  the  initial  condition  is  X-,^(o)  =  0,  then  c^  =  0  so  that 
X.,(t)  =  -^.^.(^.ate^''^i^%*).  (51) 


^12 


2:   •^5V) 


Equation  (51)  constitutes  the  second  approximation.  Thus 
far,  then 


\{t)   =  X3^3^(t)  +  X^^it) 

2  2 
=  .3(1  -  e  -*  )  -  -i^^(l-  .  2te  -*  -  i^  ^ 


(52) 


From  Eq .  (45)  the  following  is  obtained  by  substituting 
Eq.  (49)  and  (5I) 

(ttxt  -  M.i,}X,-,(t)  =  -^-j—  !i;:(l  -  e   j 5ii —  (1  -  e   }  (- 


'^   ^4 '^15 


21 

M,ut    ,    2!i2,t 
+  2te  ^  -  i-  e  ^  ). 


M.2,. 


(55) 


Multiplying  through  by  e"^4     and  collecting  terms  Eq.    (53) 
becomes 


3  3 
^    ^  -   f^  (e     ^--3+36  ^     -  e     ^    ) 


234 

2IM-4 


-tXht  M..t      2ii|,t 

e  +  (2!22^t   -   l)-(2ni^t  +1)e      +e 


(54) 


-15- 


Integration  of  Eq.    (54)   results  in 

-u.>t  -,      U-t.  t 


"3  3 
X-L3(t)e     ^     =A^(- 


2  3  4 


+  ^e 


-Ii4t 
2n^t 


1    ^H^s 


+  jx^t    -   t   -  2m,^-   -j--^(p.^t-l) 


1^4  [i-l^ 


+    C- 


(55) 


If  the  initial  condition  X-j^^(O)  =  0,  then 


=5  =  - 


2  3  4      3  3 


Thus 


Xft)   =   -2^  (  . 


13' 


3  3 

e     ^ii^ 


31 


234^ 

2p-„ 


1^4 


1  ^   ^2  »^4^  ^/4^ 


SM-i^t 


2  ^(^^4^-  -  1)  -  ^TTT^  -^  ' 


3[X4t  1 


1^4 


(56) 


Equation  (56)  constitutes  the  third  approximation  and  the 
solution  is  now 


X;^(t)  =  X^^(t)  +  X^g^^)  +  X^^(t) 


(57) 


The  process  can  be  continued  indefinitely.  Although 
involved,  the  procedure  is  straightforward. 

Referring  now  to  Eq.  (35) j  integration  of  this  expression 
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with  respect  to  t  yields 
X 


(t)  =  X^(t)  +  K^  )   sinh  a-j^X^  dt. 


Now 


sinh  a^X^  =  a^X^^  +  ■i/  +  ^,   +  .  .  .  . 

=  "l^^ll  ^  \2   ^  ^13^  ^  51(^11^^12^13)^ 

■^  51^^11  ■*"  ^12  "^  ^13)^  +  .  .  .  .         (58) 

Suppose  that  only  the  first  two  terms  of  this  expansion  are 
retained.  Then  i. 

X(t)  =  X3^(t)+K^a3^r  (X-^^+X^2+^^^^^^  ^S.   5t|  (X^^^+X^g'^X-^^^^'^^   ^^^^ 

o  o 

The  calculation  is  again  tedious,  but  can  be  carried  out  by 
anyone  acquainted  with  the  calculus , 

Equations  48  through  57  are  given  only  as  examples  of  the 
method.  The  boundary  conditions  employed,  however,  are  not 
satisfactory  as  a  solution  of  the  problem  of  Fig.  3.  Suitable 
initial  conditions  are  X^{0)  =  X^^,  XgCo)  =  X^q  and  X^q  =  X^q  = 
the  deflection  which  would  occur  in  the  springs  as  a  result 
of  P^  if  the  dash  pots  vjere  rigid. 
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Non-llnear  Maxwell  and  Volgt  Units  in  Series  with  Tvjo  Masses 

The  model  considered  in  this  section  is  shown  in  Pig.  3  and 
is  of  the  same  form  as  that  discussed  in  the  first  section 
except  that  the  elements  of  the  model  are  non-linear.  For  the 
same  reasons  given  in  the  previous  section  the  non-linear 
behavior  of  the  springs  and  dashpots  will  be  assumed  to  be  repre- 
sented by  hyperbolic  sine  functions  as  follows: 

^1  ~  ^1  ^^^  ^T^^l  ^°^  spring  1,  (60) 

Xp  =  bo  sinh  fo'^^p  ^^^   spring  2,  (6l) 

X^  -  X^=  K-j^  sinh  f-z/B-^  for  dashpot  1,  (62) 

X^  -  X^=   Kg  sinh  f^/S^  ^or  dashpot  2,  (63) 

where  b,,  b^,  a,,  ap,  are  constants  of  the  springs  and  K-.^  K^, 
6n,  pp  are  constants  of  the  dashpots. 

By  forming  the  inverse  of  Eq.  (6O-63)  the  forces  may  be 
determined  and  used  in  forming  the  equations  of  motion  as  follows: 

a-j^  sinh  \/^i   "  Q^slnh'^^iX^-X^) /K-^   =  0  (64) 

ag  sinh""-^(X^-X2)/b2  +  P>2Slhh''-^ {Xyi^) /K^-p^.slnh"-^ {i^-X^) /K^ 

F(t)-a2Sinh"^(X^-X2)/b2  -  S2Sinh"^(X^-X^/K2  =  MgX,         (66) 

Substituting  Eq.    (64)    in   (65)   and    (65)    in   (66)    the 
following  results: 

P(t)   =  m^X^  +  m^Xg   +  a^  sinh'-^X^/b^  (67) 
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from  which   X-^  =  b^  slnh  (F(t)  -  m^X,  -  m^Xg)^^  (68) 

Let  rrip  =  iD^  =0,  which  amounts  to  assuming  that  the 

accelerations  are  small  enough  to  be  neglected.   Then  X-,  is 

known . 

X-,^  =  b^  sinh  F(t)/a-^  (69) 


Also 


b-, 

Xn  =  F(t)-i  cosh  F(t)/a 
1      a.^ 


(70) 


F 


Substituting  Eq.  (65)  in  (66): 
(t)  =  mgX^  +  m^Xg  +  6^  sinh''-'-(X2  -  ^1)/% 


(71) 


Then  substituting  X,  from  Eq.  (66)  and  setting  mp  =  m,  =  0 


F(t)  =  B^sinh  "^(Xg  -  ^(t)^^  ^osh  F(t)/a3_)/K^ 


(72) 


or 


L  =  K,  sinh  ^ii  +  F(t)-i  cosh  ^^ 
t-  X  o-«         a-j       a-i 


Then  Xp  can  be  determined  as  a  function  of  F(t) 

t  b   /""^ 

X^  -  K.  1  sinh  ^i^  ^t  +  a^  / F(t)  cosh  ^^  dt 


^1 


'o 


(75) 


Substituting  Xo  and  X     in  Eq.    (66): 


2 

.tr 


F(t)   =  a^   sinh"^(X,   -  T  [k,    sinh  ^^  +  F(t)-i  cosh  ^^  dt)/b. 
^  J?      y      j_   J.  p-j^  a-j^  a^ 

o 
■^(X^   -  K^^sinh  ^  -  F(t)^  cosh4^)/K2    (7^) 


+  60   sinh 


• 

For  creep  F(t)   =  F^,   a  constant  for  t  >    0.      Thus  P(t)   =  0. 
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Then  Eq.    (7^)  becomes 


-1 


F     =  ao  slnh'-^CX^-tK.    sinh  F^/e>^)/h^ 
o  2  .       ^    p        1  oil 


-1/,; 


+  ^2   sinh"-^(X^   -  K^sinh  F^/d^)/K^ 


(75) 


The  solution  of  this  differential  equation  will  yield  the  creep 

deformation  X-,  as  a  function  of  time  t  at  a  constant  force  P  . 
5  o 


X^  =  0. 


For  relaxation  X^  =  X  ,  a  constant,  for  t  >  0.  Thus 

^  o 

Then  Eq.    (7^)   becoires 
F(t)   =  a^   sinh'YXo   "  f     %   sinh  ?^  +  F(t)^osh^lt^'  <it]/b^ 

1  J.  1   -' 


The  solution  of  this  equation  will  yield  the  force  F(t)  as  a 
function  of  time  t  at  a  constant  deformation  X  of  the  model. 
An  approximate  solution  of  Eq .  (75)  ^nay  be  obtained  as 
follows : 
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From  equation  (75) 

-  tK-^  sinh  ^q/^i)/^i  (77) 
or 

(X^  -  K-^  sinh  F^/^^)/K^   =  sinh(P^/P2 

-  a2/p2  sinh'-^(X^  -  tK-^  sinh  Pq/^i)/^i  (78) 

Expanding  the  hyperbolic  sine  in  the  right  hand  side  of  (78)  and 
substituting  in  (78) 

X^  =  K^a^  +  K^cxg  ^Qg^  j^^/^^  sinh'-'-CX^-  tK-La-j^)^^^! 

-  KgQg   sinhja^/Sg   sinh'-^^CX^   -  tK^a^)/b-j_  1  (79) 

where  a^    =  sinh  Pq/Ptj     o'p  ~  sinh  Pq/^o  ^^^  *^3  ~  cosh  F  /pp. 

Let    (X^  -  tK-,a^)/b-,    =  Z-^;    then  X^  =  h-,Z^  +  K-,a-, 
and   (79)  becomes 

Z,   =  Kgttg/b^  cosh    \^2l^^2.   ^^^^"  ^^ 

-  K^a^/b-L  sinhf 32/62   sinh""*"  Z,  1  .  (80) 

Using  the  series  expansion  of  the  inverse  hyperbolic  functions 
and  dropping  terms  above  the  fourth  order  equation  (80)  is  found 
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to  be  -2 


(81) 

Also  using  the  series  expansion  of  the  hyperbolic  functions: 

•7. 


coshja^/B^lZ^  -  J-)j- 


1  +  -^r~"  ^3  ^  ^"Ti 5r"^^3  ^"72 W^^l> 

4/4      4/4       4/4 
^   ^2^^2  „8   ^2^^2  „10  ^  ^2/^2  „12  .,^  , 

■*■   144   ^3  "  129b   ^-^   ^  ynm  ^3  •  182) 


and      _  7?  -x  Z^ 


3-o/^o    ^    Z^    Z-r      Z 


,-2^(z|.^.3|-  ^)  (83) 

Substituting  (82)  and  (83)  in  equation  (8I),  dropping 
terms  above  the  fourth  order,  and  transposing,  (81)  becomes 

Z^  +  Yi^^/\^^{^^/^^)Z^   -  K2a2/b^(-|-f-^)z| 

a|/3^    a  /Pp  , 
+  K^a^/b^C-^^-j^ 2^)Z^ 

^\/^\        a|/3^  ,  (^^) 

Equation  (84)  is  now  in  a  form  to  which  the  method  of 

■5 
reversion^  can  be  applied.   To  simplify  the  computations  let 

p   ^      4 
&l»  &2'   ^3  ^^^  ^4  ^®  ^^^  coefficients  of  Z,,  Z^, Z  ^  and  Z^ 
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respectively,  and  let  KgOtp/t*-!  =  B. 
Then  Eq.  (84)  becomes 

Z^  +  g^Z^  +   ggZ^  +  g^Z^  +  gi^Z^  =  B.  (85) 

This  expression  may  be  solved  in  the  same  manner  outlined  on 
page  12  for  the  solution  of  Eq.  42. 

The  boundary  conditions  which  may  be  applied  to  solve  this 
problem  for  creep,  i.e.  F(t)  =  F  =  constant,  are: 

X^(0)  =  X2(0)  =  X,(0)  =  b^  sinh  'Pq/s-^   and 


X^{0)   =  X^iO)   =  0.  (86) 

For  stress  relaxation,  i.e.,  X^(t)  =  X,q  =  constant,  the 
boundary  conditions  which  may  be  employed  are: 

X  (0)  =  X^{0)   =  X^Q,   F(0)  =  a^sinh"-'-  X^^/b^^,  and  ^2(0)  =  0. 
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SUMMARY 

The  method  of  predicting  creep  in  bending  from  data 
on  creep  in  tension  previously  described  has  been  extended 
to  the  prediction  of  creep  in  bending  when  the  creep  in  ten- 
sion and  compression  are  iinequal  and  when  the  time -dependent 
and  time -independent  stress  functions  are  unequal.   It  was 
shown  that  the  stress  distribution  and  position  of  the  neutral 
axis  in  a  beam  changed  with  time  when  the  creep  in  tension 
and  compression  were  unequal;  and  the  stress  distribution 
changed  with  time  when  the  time -dependent  and  time -independent 
stress  functions  were  unequal. 

Computed  creep  deflections  were  compared  with  avail- 
able data  on  canvas  laminate  and  polystyrene,  with  good 
results  for  the  former  and  fair  results  for  the  latter. 


CONCLUSIONS 
It  has  been  shown  that  creep  and  stress  distribu- 
tion in  bending  can  be  predicted  from  tension  and  compres- 
sion creep  data.  However,  the  available  test  data  are 
inadequate  to  provide  an  accurate  check  of  the  theory. 


ACKNOWLEDGEMENTS 

This  project  was  conducted  in  the  Department  of 
Theoretical  and  Applied  Mechanics  as  part  of  the  work  of 
the  Engineering  Experiment  Station  of  the  University  of 
Illinois  in  cooperation  with  Picatinny  Arsenal,  Ordnance 
Corps,  Department  of  the  Army. 

The  authors  are  grateful  to  B.  K.  Ghandhi, 
W.  A.  Hagemeyer,  G.  Khosla  and  H.  G.  Russel  for  assis- 
tance in  the  computations  and  preparation  of  illustrations 


INTRODUCTION 

When  a  material  is  subjected  to  a  constant  stress  the  time- 
dependent  portion  of  the  resulting  strain  imposed  on  it  is  referred 
to  as  creep.  Although  creep  behavior  of  materials  under  vari- 
ous states  of  stress  is  both  of  practical  and  theoretical 
interest,  most  creep  data  reported  in  the  literature  have  been 
obtained  with  simple  tensile  loading.   In  the  design  of  many 
industrial  products  one  of  the  more  simple  cases  of  complex 
stresses  encountered  is  bending.  The  purpose  of  the  present 
investigation  was  to  develop  a  theory  correlating  the  creep 
behavior  in  bending  with  creep  in  tension  and  compression. 

In  a  previous  paper  (l)*  the  existing  literature  re- 
lated to  this  problem  was  reviewed  and  a  method  for  predicting 
creep  in  bending  from  the  creep  data  in  tension  was  derived. 
In  the  previous  paper  the  creep  behavior  of  the  material  was 
assumed  to  be  equal  in  tension  and  compression  and  the 
coefficients  of  the  time -dependent  and  time -independent  terms 
of  the  stress  function  were  considered  equal.  The  theoreti- 
cal equations  obtained  from  this  analysis  were  compared  with 
experimental  data  obtained  from  a  canvas  laminate  with  good 
success.  The  method  was  also  applied  to  the  case  of  creep 
under  a  nonuniform  bending  moment. 

An  earlier  paper  by  Marin  (2)  which  has  not  been  re- 
ferred to  in  recent  papers  was  discovered  since  preparation 

*  Numbers  in  parentheses  refer  to  the  list  of  references 
appended  to  this  paper. 
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of  our  first  paper  on  this  subject.   In  the  paper  by  Marin 
creep  was  assumed  to  be  described  by  a  power  f\inction  of 
time  and  the  stress  functions  in  the  time -dependent  and 
time -independent  terms  were  also  assumed  to  be  power  func- 
tions. The  predictions  of  Marin *s  analysis  were  tested 
against  creep  data  from  four  aluminum  alloy  beams  with  errors 
of  13,  22,  24  and  28  per  cent,  respectively.  The  application 
of  the  analysis  to  statically  indeterminate  structures  was 
also  described. 

The  present  paper  extends  the  analysis  to  include 
bending  of  materials  which  have  unequal  creep  in  tension  and 
compression  and  which  exhibit  unequal  coefficients  of  the 
time -dependent  and  time -independent  terms.  Both  of  these  con- 
ditions make  the  solution  of  the  problem  much  more  difficult. 

BENDING  WITH  UNEQUAL  CREEP  IN  TENSION  AND  COMPRESSION 

Derivation: 

Some  materials  evidently  have  imequal  creep  behavior 
in  tension  and  compression  (5).  The  derivation  of  creep  re- 
lations for  bending  from  unequal  relations  for  tension  and 
compression  creep  is  complicated  by  the  lack  of  symmetry  which 
probably  causes  the  stress  distribution  and  position  of  the 
neutral  axis  to  change  with  time. 

In  previous  papers  (4,  5,  6)  the  creep  curves  ob- 
tained from  a  canvas  laminate  tested  at  different  constant 
stresses  in  tension  were  found  to  be  closely  described  by  an 
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equation  of  the  form 

€  =  €  +  mt'^  =  (e  '  +  m't'^)  sinh  q/q^  (l) 

o  o  ^ 

where  €  is  the  strain;  t  is  the  time  under  a  constant  stress, 
a;  €*,  m'j  n  and  a  are  constants  whose  values  depend  on  the 
material  and  temperature  ^  and  e  and  m  are  functions  of 
stress. 

The  stress  dependence  in  this  relation  is  of  the  same 
form  as  that  in  the  activation  energy  theory  of  creep  as 
advanced  by  Kauzmann  (7).   It  has  also  been  shown  (6)  that 
an  equation  for  the  stress -strain  curve  obtained  from  a  con- 
stant strain  rate  tension  test  could  be  derived  from  equa- 
tion (1)  and  that  the  resulting  equation  was  in  close  agree- 
ment with  the  experimental  data.  It  should  be  noted  that 
the  time-dependent  term  €  =  €'  sinh  cj/a  in  equation  (1) 
expresses  the  instantaneous  strain,  i.e.,  the  strain  at 
time  zero,  and  that  it  does  not  represent  linear  elasticity 
except  when  a^  is  large  compared  to  a. 

Equation  (l)  accurately  represents  the  creep  data 
of  some  materials  under  conditions  of  constant  stress.   It 
probably  does  not,  however,  accurately  describe  creep  under 
varying  stress.  Thus,  if  the  stress  distribution  in  the 
beam  changes  during  creep,  as  described  by  Popov  (8)  and 
Pried  (9)  some  inaccuracy  may  be  introduced  by  the  use  of  this 
equation  unless  the  material  obeys  the  time -hardening  law  (lO) 
which  is  probably  not  the  case . 
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If  the  constants  in  the  equation  for  tension  are  all 
different  from  those  for  compression  the  creep  equations  may 
be  written 

Crp  =  (€q;j,  4-  m'^t  )sinh  (s.j^/o^,j,  (2) 

S  =  (^is  -^"^'s^''^)^^"^  ^s/^os  (5) 

where  the  subscripts  T  and  s  denote  tension  and  compression 
respectively. 

To  determine  the  bending  creep  relation  the  conditions 
of  equilibrium  and  geometry  must  be  satisfied.  From  the 
observation  that  plane  sections  remain  plane  in  a  beam  sub- 
jected to  a  uniformly  or  nearly  uniformly  distributed  bending 
moment  the  following  well-known  geometrical  relations  are 
obtained : 

s   p  *   T   p  ^  ' 

where  yrn  and  y  are  the  normal  distances  from  the  neutral 
axis  to  any  point  in  the  tension  or  compression  side  of  the 
beam,  respectively,  €„  and  e  are  the  corresponding  strains 
and  p  is  the  radius  of  curvature  of  the  neutral  axis .  The 
neutral  axis  is  at  some  undetermined  distance  c  from  the 
extreme  "fiber"  on  the  tension  face  of  the  beam. 

The  fact  that  plane  sections  in  a  beam  subjected  to 
pure  bending  remain  plane  has  been  demonstrated  for  plastically 
bent  beams  by  Nadai  (11)  and  for  creep  of  lead  beams  by 
MacCullough  (12).  Also,  consideration  of  compatibility  of 
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strains  in  the  different  "fibers"  of  a  beam  of  great  length, 
bent  by  a  constant  bending  moment  into  the  arc  of  a  circle, 
suggests  that  plane  sections  must  remain  plane  at  all  points 
along  the  beam  which  are  remote  from  the  ends.   If  this  were 
not  so, different  "fibers"  would  be  markedly  displaced  rela- 
tive to  each  other  at  points  at  a  distance  from  the  center  of 
the  beam. 

The  conditions  of  equilibrium  for  a  beam  are  expressed 
by  the  following  equations  when  the  accelerations  involved 
in  the  creep  are  so  small  as  to  be  negligible — which  is 
usually  the  case: 

/   adA  =  0  (5) 

.c 
aydA  =  M  (6) 

c-h 

where  A  is  the  cross -sectional  area  of  the  beam,  h  is  the 
depth  of  the  beam  and  c  is  the  distance  from  the  neutral 
axis  to  the  tension  face  of  the  beam. 

Solving  for  a^  and  a^   in  equations  (2)  and  (5)  and 
substituting  equations  (4)  for  e  and  €«  the  following  equa- 
tions are  obtained. 

-1   yrn/f  n  N^ 

0^   =  a^^sinh  —   =  a^^sinh-  ^  y^   (7) 

^oT  "^  "^*T^ 
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^s  =  °os  slnh-l— IS-^ =  c^3  sinh-1  f  y^  (8) 


Where  the  symbols 


N  =  1 __  ;  and  N„  = 


T  n^  '  "-    ^'s     ,       rig 

have  been  introduced. 

For  a  beam  of  rectangular  cross -section  of  width  b, 
dA  =  bdy.  Making  this  subsitution  equations  (5)  and  (6) 
become:  q  c 


■bj       a^pdy  +  b 
-(h-c)       "0 


bJ   a^c 


(Jgdy  =  0  (9) 


and 


M  =  b/    a^y^dy  +  b/  cT^y^dy       ,     (10) 
-(h-c)        ^0 

Substituting  for  a^  and  a_  from  equations  (7)  and  (8) 

x  S 

into  equation  (9),  and  integrating,  the  resulting  expression 
becomes: 


^[x4(.-c)s.n.-;i(.-c)./77^^^^^     ^^ 


jj2i[l.^csl„h-l^c-Jl.(p^o)2(.  (11) 


Substituting  for  a^n  and  a^  from  equations  (7)  and  (8) 
into  equation  (10)  and  integrating 
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M  = 


2  Cr 


N„ 


(h  -  c)^  +  1 


-1  ^T, 
sinh    ~(h  -  c) 


-^(h-c)  Jl  + 


Nn 


•(h-c) 


-1  ^s 
c[  +  U  sinh   ~  c 


(12) 


For  a  given  constant  bending  moment  M  applied  for  a  given  time 
t  to  a  given  beam,  equations  (11 )  and  (12)  are  simultaneous 
equations  in  two  unknowns,  the  location  of  the  neutral  axis 
c  and  the  radius  of  curvature  p .  Values  of  p  and  c  may  be 
obtained  by  trial . 

For  a  beam  of  constant  bending  moment  the  deflection 
at  mid-span  may  then  be  obtained  from  the  well  known 
geometrical  relation: 


z  = 


5p 


(15) 


Where  Z  is  the  deflection  at  the  center  relative  to  the  ends 
of  a  span  of  a  beam  of  length  H   subjected  to  a  uniform  bending 
moment . 

The  stress  distribution  may  be  obtained  from  the  value 
of  c  and  equations  (7)  and  (8). 

Since  N  and  Nm  are  functions  of  time  the  radius  of 
curvature  p  will  change  with  time.  Since  Ng  and  N^  are 

unequal  functions  of  time  it  is  evident  that  the  ratios 

^s      N^ 

-—  and  —  cannot  both  remain  constant  and  independent  of  time 

r  r 

Hence,  from  equation  (ll)  it  follows  that  c  must  vary  with 
time. 
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Since  the  ratios  —  and  — =-  change  with  time  the 

P      P 

stress  distribution  must  change  with  time.  This  follows  from 
equations  (7)  and  (8).  Not  only  does  the  stress  distribu- 
tion change  with  time  because  of  the  time  dependence  of  the 
above  ratios  but  because  the  position  of  the  neutral  axis 
changes  with  time  and  this  affects  the  values  of  y™  and  y 
in  equations  (7)  and  (8). 

Equation  (l)  may  accurately  describe  the  relation- 
ship between  strain,  stress  and  time  only  for  constant 
stress  or  stresses  which  change  only  slightly.  Hence,  if 
the  change  in  stress  distribution  described  above  causes 
large  changes  in  stresses  during  creep  then  the  above  deriva- 
tion will  not  yield  accurate  results . 

It  is  easily  verified  that  equations  (ll)  and  (12) 
reduce  to  equations  (l6)  and  (18)  in  the  previous  paper  (l) 
when  tension  and  compression  creep  are  equal.  This  verifica- 
tion is  accomplished  by  substituting  the  following  in  equations 
(11)  and  (12): 

^oT  =  ^os  =  ^o'   N^  =  Ng  =  Np;  h  -  c  =  c. 

Illustration: 

Equations  (ll)  and  (12)  have  been  solved  for  the  par- 
ticular problem  of  a  beam  of  canvas  laminate  in  which  the 
creep  in  tension  and  compression  can  be  described  by  relations 
of  the  form  of  equations  (2)  and  (3).  And  the  constants  for 
creep  in  tension  were  the  same  as  employed  for  the  canvas 
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laminate  analyzed  in  the  previous  report  (l); 

^oT  "  ^°°^  ^^^ 
n^  =  0.1185 

while  the  constants  for  creep  in  compression  were  arbitrarily 
chosen  to  be  20  per  cent  higher  than  those  for  tension  creep; 

€   =  m •  =  0.002250 
OS     s 

a^^  =  4800  psi 
OS         ^ 

ng  =  0.14196 

To  solve  equations  (11)  and  (12)  for  the  purpose  of 
determining  the  deflection-time  relation  in  bending  the 
values  of  N  and  Tm  were  computed  for  a  particuleir  set  of 
values  of  time  t.  Then  the  constant  coefficients  in  the  equa- 
tions were  computed.  Graphical  techniques  for  solving  the 
simultaneous  equations  were  tried  but  were  found  to  lack  the 
necessary  precision. 

A  numerical  procedure  was  employed  in  which  Nm(h-c)/p, 
and  N  c/p  were  used  as  variables.  Likely  values  of  these 
variables  were  inserted  in  equation  (ll)  and  the  value  of  the 
left  hand  side  and  the  right  hand  side  of  the  equation  were 
computed.  Then  the  values  of  both  variables  were  plotted  as 
ordinates  versus  the  value  of  the  right-or  left  hand  side  of 
the  equation.  Thus  any  number  of  pairs  of  values  which  would 
satisfy  equation  (ll)  could  be  read  from  these  curves.  Pairs 
of  these  values  were  then  substituted  in  equation  (12)  and 
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the  bending  moment  was  computed.  This  process  was  repeated 
imtil  the  computed  bending  moment  matched  the  required  value 
with  sufficient  accuracy. 

From  the  values  of  the  above  variables  which  satis- 
fied equations  (ll)  and  (12)  the  values  of  p  and  c  were 
computed,  and  from  p  the  deflection  was  computed  by  means 
of  equation  (13).  The  deflection-time  relation  for  a  bend- 
ing moment  of  679  inch-pounds  determined  in  this  way  is 
shown  in  Pig.  1  together  with  the  corresponding  experimental 
bending  creep  curve  and  the  theoretical  creep  curve  for 
equal  creep  in  tension  and  compression. 

The  follovjing  observation  was  of  interest:   In  spite 
of  the  fact  that  the  assumed  creep  in  compression  varied  from 
20  per  cent  greater  than  in  tension  at  t  =  0  to  35  per  cent 
greater  than  in  tension  at  880  hours  for  the  same  stress,  the 
computed  creep  in  bending  was  nearly  the  same  over  the  time 
interval  considered  as  that  computed  for  a  beam  in  which  the 
creep  in  compression  was  the  same  as  that  in  tension. 

The  position  of  the  neutral  axis  was  found  to  vary  with 
time  from  one  side  of  the  center  line  to  the  other  as  shown 
in  Pig.  2. 

The  distribution  of  stress  across  the  beam  was  com- 
puted at  two   different  values  of  time  from  equations  (7)  and 
(8),  and  was  plotted  for  a  time  of  zero  and  of  880  hr.  in 
Pig.  3a.   It  was  observed  that  the  change  in  stress  distribu- 
tion was  relatively  small — less  than  2  per  cent  in  880  hours 
at  the  extreme  fiber. 
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In  order  to  better  illustrate  the  change  which 
occurred  in  the  stress  distribution  the  difference  between 
the  computed  stress  distribution  and  the  equivalent  linear 
stress  distribution  was  determined  for  three  values  of  time 
and  is  shown  in  Pig.  5b.   It  was  found  that  a  high  degree 
of  precision  was  required  in  the  computations  to  insure 
accurate  results  in  the  stress  distribution* 

BENDING  WITH  UNEQUAL  TIME  DEPENDENT  AND  INDEPENDENT 

STRESS  FUNCTIONS 

Derivation: 

The  creep  behavior  of  some  materials  is  accurately 

represented  by  a  more  general  form  of  equation  (l)  in  which 

the  constants  in  the  stress  function  are  different  for  the 

time -independent  and  time -dependent  terms  as  follows: 

€  =  e'  sinh  q/o^   +  m't^sinh  a/a_  (l4) 

o        e  m  ^      ' 

When  the  tension  and  compression  creep  behaviors  are  identi- 
cal and  are  single -valued  functions  of  stress  and  time  as 
described  by  equation  (l4)  for  constant  stresses  or  small 
changes  in  stress,  the  creep  in  bending  can  be  expressed  by 

equations  derived  in  the  following  manner: 

Equation  (l4)  may  be  written  as  follows  by  dividing  by  e^ 
e  =  sinh  x  +  T  sinh  ax  (15) 

where 

,^n 

e  =  €/€q  ;  T  =  ^^-^   ',     a  =  a^/a^  and  x  =  <5/(5^.        (l6) 

For  a  straight  beam  of  rectangular  cross  section  in 
which  the  creep  equations  in  tension  and  compression  are 


....-nein^  ,fS?  ffl  \  'iWf''' 
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identical  the  neutral  axis  can  "be  shown  to  be  at  half  the 

depth  of  the  beam.  Thus  from  statics  the  bending  moment 

M  in  a  beam  of  width  b  and  depth  h  is  related  to  the  stress 

a   by  the  following  equation: 

c  c 

M  =J   ho-^jdj   =  by   (Jg  xy  dy  (17) 

-c         -c 

where 

c  =  half  the  beam  depth,  and 

y  =  distance  measured  from  the  neutral  axis  in  a 
direction  perpendicular  to  the  neutral  axis . 

Prom  the  geometry  of  strains  as  expressed  by  equation 
(K)   and  the  creep  relation,  equation  (l5)j  the  following  was 
obtained: 

y  =  pe  =  p6  fe  =  pe'(sinh  x  +  T  sinh  ax).    (l8) 

Differentiating  equation  (l8)  at  constant  p  and  T 

dy  =  pe  '  (cosh  x  dx  +  aT  cosh  ax  dx)       (19) 

Substituting  equations  (l8)  and  (19)  into  equation 

(17)  for  any  given  time  t,  (and  corresponding  T  and  p)  and 

expanding  r   x-, 

M  =  ba  (pe  ')  <  /   x  sinh  x  cosh  x  dx 


-^1 

+  aT  /   X  sinh  x  cosh  ax  dx  +  T/  xs3nhax  coshx  dx 


/-I 
V 


2^' 

+  aT  y   X  sinh  ax  cosh  ax  dxj  (20) 

-^1 
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where  the  limits  of  x  correspond  to  the  limits  of  y;  that 

is,  Xn  is  the  stress  at  y  =  c  and  -x-,    is  the  stress  at  y  =  -c. 

Performing  the  integration  as  indicated  in  equation 
(20),  the  resulting  equation  for  radius  of  curvature  p  as  a 
function  of  bending  moment  M  and  stress  at  the  extreme 


"fiber"  a^  =  ^^^^i   was 


ba^(e„'  )2p2. 


M  =s 


(2  sinh  x^  +  l)x^  -  sinh  x^  cosh  x^ 


+  ^1(2  sinh  ax,  +  l)ax-j^  -  sinh  ax-j_  cosh  ax^ 

+  ^TJ  x.  sinh  x^    sinh  ax, -^ (a  sinh  x, cosh  ax. 

L  1       ^    ^    ^    (a2  -  1)  ^      ^ 


-  cosh  X-,  sinh  ax,  ) 


(21) 


Solving  equation  (l8)  for  l/p  when  y  =  c  and  x  =  x^  the 
following  is  obtained: 

±   =  _^(sinh  X,  +  T  sinh  ax,  ).         (22) 
p    c  ^      1  1' 

These  two  independent  relations,  (21)  and  (22),  can  be 
solved  simultaneously  at  any  given  time  t  =  t,  (or  T  =  T, ) 
for  the  extreme  fiber  stress  function  x-,  and  the  radius  of 
curvature  p.  The  deflection  for  a  beam  having  uniform  bending 
moment  may  then  be  obtained  from  equation  (15). 

If  the  ratio  a  -  ^^^^^     is  allowed  to  equal  unity, 
equation  (l4)  reduces  to  equation  (l).   Consequently, 
equation  (21)  should  reduce  to  equation  (l8)  of  reference  (l). 
Performing  the  indicated  substitution  and  solving  for  M 
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M 


ba  (p€  ')^f      o  r      P 
=  — P   9 — }(i  +  T  )  (2  sinh  x^  +  1)  x-j^  -  sinh  x-j^cosh  x^ 


+  T 


x.  sinh^x.  -  -^^   {-o-^ — )(a  sinh  x,  cosh  ax, 
^      ^       a — >1  a"^  -.1         ^ 


-   cosh  x^  sinh  ax^ ) 


.(25) 


Taking  the  limit  and  making  use  of  equation  (l4)  of  reference 
(1)  equation  (23)  V7as  found  to  reduce  to  equation  (l8)  of  refer- 
ence (l). 

The  simultaneous  equations  (21)  and  (22)  may  be  solved 
by  substituting  equation  (22)  for  p  in  equation  (21),  expand- 
ing and  collecting  terms  in  T.  The  following  quadratic  equa- 
tion in  T  results: 

T^  j(2MA>c^a^)sinh^ax-j^  -  x^(l  +  2  sinh^ax-j^) 
+  (l/a)  cosh  ax^sinh  3X3^+  tI   (4M/bc^-  kx^ 

4/(a^-  l)Jsinh  ax^cosh  x^  +  Ua/(a^-  l) 
+  1  sinh  x-|^cosh  x-^    +  (2M/bc  a  )sinh  x^  -  Xn(l  +  2  sinh  x-^  =  0 

(24) 
Equation  (24)  is  then  evaluated  for  the  given  value  of  M  and 
several  reasonable  values  of  x.  by  using  the  quadratic  formula. 
The  solution  of  the  quadratic  formula  gives  the  time  function 
T  corresponding  to  each  chosen  stress  function  x^.   From  the 
value  of  T  the  time  t  was  determined  from  equation  (I6).  The 
curvature  l/p  was  determined  from  equation  (22)  by  substituting 
T  and  x-j^.  Then  by  utilizing  equation  (13)  the  deflection  z  was 
was  determined  for  the  same  value  of  time. 


sinh  ax-,  sinh  x. 


cosh  ax^slnh  x. 
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The  above  procedure  is  straight  forward  "but  involves 
a  great  deal  of  labor  and  requires  a  large  number  of  signifi- 
cant figures  in  the  computations. 

From  the  values  of  x-,  at  different  times  the  change 
in  extreme  fiber  stress  may  be  determined  as  a  function  of 
time  by  use  of  the  last  relation  in  equations  0-6).  Also 
at  any  given  value  of  time  and  corresponding  values  of  T  and 
p  the  stress  distribution  can  be  computed  from  equation  (l8). 
Because  the  creep  functions  were  assumed  to  be  the  same  in 
tension  and  compression  the  neutral  axis  does  not  change 
position  with  time. 

Illustration: 

Creep  data  were  available  for  tension  {13 >   1^)  and 
bending  (13)  of  polystyrene.   The  two  sets  of  tension  creep 
data  were  both  made  on  polystyrene  of  the  same  description 
from  the  same  manufacturer.  One  set  of  tension  creep  data 
did  not  have  enough  self  consistency  to  permit  evaluating 
the  constants  of  the  creep  equation,  but  by  combining  the  two 
sets  of  data  this  became  possible — though  with  limited 
accuracy. 

It  vjas  found  that  equation  (l)  could  not  adequately 

describe  the  data  so  equation  (l4)  was  tried  with  the  results 

shown  in  Fig.  4.  The  corresponding  values  of  the  constants 

in  equation  (l4)  were  found  to  be  as  follows: 

^Q   =  4. 3077  X  10"^ 

a^   =  20,000  psi 

m'  =  4.098  X  10"^ 

o^   =  700  psi 
m       ^ 

n  =  0.475 
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These  values  were  obtained  largely  by  trial.  The 
data  did  not  yield  readily  to  rational  procedures  of  evalu- 
ating the  constants  such  as  used  previously  (4,  6).  Utiliz- 
ing these  values  of  the  constants  and  the  procedure  outlined 
above,  the  creep  in  bending  was  computed  for  a  beam  of  the 
dimensions  employed  in  the  tests  of  reference  (15)  and  a 
bending  moment  of  553  in.  lb.  The  results  of  these  computa- 
tions are  shown  in  Pig.  5  together  with  test  data  reported 
in  reference  (13). 

Also  shown  in  Pig.  5  is  the  change  in  the  stress  at 
the  extreme  fiber  with  time.  The  computed  stress  at  the 
extreme  fiber  decreased  10  per  cent  in  1000  hours  and  de- 
creased 12.4  per  cent  in  3200  hours.  Since  the  bending  mo- 
ment remained  constant  the  stress  distribution  changed. 

The  agreement  between  the  theoretical  and  experimental 
creep  deflections  was  not  as  good  as  would  be  desired 
although  the  shape  of  the  curve  was  satisfactorily  predicted 
from  100  to  1000  hr.  Possible  reasons  for  the  lack  of 
better  agreement  may  be  among  the  following:  (a)  Creep  of  the 
polystyrene  tested  in  bending  and  part  of  the  tension  tests 
may  not  have  been  accurately  described  by  the  constants 
selected  for  the  creep  equation.   Some  of  the  tension  creep 
curves  for  this  material  shov/ed  similar  departures  from  the 
theory  in  both  tension  and  bending  tests;  (b)  crazing  was 
observed  to  occur  in  tension  creep  and  on  the  tension  side 
of  the  beams;  (c)  the  data  in  the  first  100  hours  of  creep 
in  bending  appear  to  be  questionable;  (d)  creep  in  tension 
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and  compression  may  not  have  been  the  same;  and  (e)  the  change 
in  stress  distribution  with  time  caused  the  stress  on  a 
given  fiber  to  change  during  creep  instead  of  remaining  con- 
stant as  assumed  by  the  use  of  equation  (l^). 
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o      Test  data,  from  reference(l) 
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— A — Theory  for  unequal  creep  in  tension  and 
compression ,  from  equations  (11,12,13) 
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THEORIES,  FOR    CANVAS    LAMINATE    AT    11^?.    AND    BENDING 
MOMENT   OF   679  IN.-LB. 
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FIG.2     LOCATION    OF    NEUTRAL     AXIS    DURING    CREEP    IN     BENDING 

WHEN     TENSION     AND     COMPRESSION     CREEP     ARE      UNEQUAL, 
FOR      CANVAS     LAMINATE      AT     77*  F     AND      BENDING       MOMENT 
OF      679  IN.-LB. 


Theory  from   Equation  (14) 

oo  i^t.  Data   from  Reference  (13) 

•  •  •  •  Data  from  Reference  (14) 

Bi  Indicates   Fracture 
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FIG.  4      CREEP   TESTS    IN    TENSION    AT   DIFFERENT 
STRESSES    FOR    POLYSTYRENE    AT  77 F 
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SUMMARY 


Results  of  creep  tests  of  the  following  plastics  are 
reported:   (a)  polyethylene 

(b)  Polychlorotrif luoroethylene  (Kel-F),  crystal- 
line, and 

(c)  polyvinyl  chloride  (Geon  ^04),  annealed. 
Each  of  the  plastics  was  loaded  in  tension  creep  at 

several  different  values  of  stress.  Creep  tests  were  performed 
at  a  temperature  of  77°F  and  a  relative  humidity  of  50^. 

Creep  strains  observed  in  these  tests  were  compared  with 
those  obtained  theoretically  by  assuming  a  creep  equation  of 
the  form: 

€  =  €'  sinh  ~  +  m't^sinh  ~ 

where  e  =  strain,  a  =  stress,  t  =  time  and  n,  e',  a  ,  m' 
and  a^  are  constants . 

It  is  shown  that  the  experimental  creep  curves  are  very 
well  described  by  the  above  creep  equation. 


CONCLUSIONS 

It  was  observed  that  a  creep  equation  of  the  i'orm 
6  =  €'  sinh  5—  +  "''■t  sinh  — - 

where  e',  a  ,  m',  a  and  n  are  constants,  described  the 
o'   €       m 

experimental  creep  data  for  the  three  plastics  very  well. 

The  rate  of  creep  at  500  hours  and  the  percentage 
increase  in  strain  from  20  seconds  to  500  hours  was  the 
greatest  for  Polyethylene  and  the  least  for  Kel-F 
(crystalline).  Hov^ever,  the  total  creep  at  500  hours  was 
the  greatest  for  Kel-F  (crystalline)  and  the  least  for 
Geon-4o^  (annealed). 

Crazing  of  Kel-P  (crystalline)  occurred  at  the 
highest  stress  of  a  =  5600  psi.,  so  that  the  total  creep, 
creep  rate  and  percentage  increase  in  creep  were  unusu- 
ally high  at  this  stress.  Further  study  in  this  connec- 
tion is  warranted. 
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AN  EQUATION  FOR  TENSION  CREEP  OF  THREE 
UNFILLED  THERMOPLASTICS 


Introduction 


Efforts  to  express  creep  behavior  by  a  mathematical 
expression  have  been  attempted  by  various  authors.  Creep 
strain  at  constant  load  has  been  described  by  both  linear 
and  non-linear  functions  of  time. 

For  the  plastics  under  consideration,  creep  strain 
is  a  non-linear  function  of  time,  as,  at  no  stage  of  the  creep 
curve  was  the  creep  rate  constant.   A  few  of  the  investiga- 
tions in  which  creep  has  been  observed  to  be  a  non-linear 
function  of  time  are  given  below: 

Leaderman^  '   expressed  creep  of  Bakelite  under  torsion 

as: 

€  =  A-log  t  +  B-t  +  C  (l) 

where  e  =  creep  strain,  t  =  time  and  A,  B,  and  C  are  constants 
which  are  functions  of  stress  and  of  the  material. 

The  same  type  of  equation  has  been  used  by  Telfair, 

(2) 
Carswell  and  Nason^  ''  for  phenolic  plastics  and  later  by 

(^) 

Lyons ^-^''  for  tire  cords. 

(h.) 
Cottrel  and  Aytekin^  '  used  a  modification  of  the 

Andrade  equation  for  single  crystals: 

e  =  €q  +  Af^^^  +  Bt  (2), 

where  €  =  shear  strain  at  time  t,   ^     =   instantaneous  shear 
strain,  and  A  and  B  are  constants. 


*  Superior  numbers  in  parentheses  refer  to  the  references  at  the 
end  of  this  paper. 


2. 


Pao  and  Marin^-^'^  used  the  following  expression  for 


creep : 

€  =  A  +  B(l~e"^^)  +  Dt  (3) 

where  again  constants  A,  B,  c  and  D  are  functions  of  stress 
and  of  the  material . 

For  the  present  tests,  however,  the  following  creep  equa- 
tion^ '*'    '  was  employed: 


€  =  €_  +  mt^^  W 


■+■   mr 
o 


where  n  is  a  constant,  while  e  and  m   are  hyperbolic  sine  func- 
tions of  stress  a   and  are  expressed  by 

S  =  ^i  slnh  f- 

(5) 

m     =  m'  slnh  — — 

m 

where  e^^  m^  ,    a  and  a     are  constants  of  the  material. 

Writing  Eq.  k   in  the  final  form  by  the  help  of  Eq.  5« 


€  =  € •  sinh  I-  +  m't^inh  |-  (6) 

€  m 

Eq.  6  has  been  used  to  describe  the  creep  of  polyethylene, 

Kel-F  (crystalline)  and  Geon-4o4  (annealed)  in  this  report  and 

shows  a  good  agreement  with  the  experimental  results. 

Materials 

Polyethylene  was  obtained  in  the  form  of  a  compression 
molded  sheet  9/16"  thick  from  Plax  Corporation,  Con.   The  mold- 
ing powder  used  in  fabricating  the  sheet  of  polyethylene  was 
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Bakellte  DE  2^00,  having  a  melt  index  of  about  1.2  decigrams/ 
minute  with  an  extrapolated  number -average  molecular  weight  of 
approximately  22,000.  There  were  no  additives  or  anti-oxidants 
employed . 

The  molding  powder  was  preheated  for  25  minutes  at  a 
pressure  of  5  psi.   The  pressure  was  then  raised  to  approximately 
30  psi.  for  ^  minutes  and  then  left  in  the  press  for  about 
30  minutes  with  water  at  60°P  running  through  the  cooling 
platens.   It  was  reported  that  such  a  thick  sheet  probably  con- 
tained a  large  proportion  of  the  crystalline  phase. 

Kel-F  (polychlorotrif luoroethylene )  of  grade  300  was 

obtained  in  sheets  of  1/8  in.  thickness  from  the  M.  W.  Kellogg 

Co.  This  grade  was  unplasticized  and  represents  the  highest 

molecular  weight  in  the  form  of  molding  powder.  The  high  density 

molding  powder  had  a  compression  ratio  of  2  to  1.   The  sheets 

were  compression  molded  under  a  pressure  of  2500  to  3000  psi. 

and  at  a  temperature  of  480  to  500°F.  The  crystalline  form  of 
Kel-F  was  obtained  after  pressing  by  dropping  the  temoerature 
quickly  by  means  of  steam  to  450^P  then  cooling  to  room  tempera- 
tMre  in  the  press.   It  required  6  to  8  hours  to  cool. 

Geon  4o4  (polyvinyl  chloride)  was  obtained  as  a  sheet 

9/16  inch  thick  from  the  Nixon  Nitration  Works .   The  formulation 

was  designated  VX-I08  Natural  and  contained  no  plasticizer, 

3.3^  dibutal-tin-diborate  as  a  stabilizer,  and  0.5^  stearic  acid 

as  a  lubricant.  The  ingredients  were  mixed  and  blended  at 

280  P  and  calendered  into  sheets  0.02  in.  thick  on  rolls  having 

a  temperature  of  300°P.  About  28  of  these  sheets  were  laminated 


together  in  a  hydraulic  press  to  form  the  final  sheet  and 
remove  most  of  the  residual  strain  resulting  from  the  calender- 
ing. The  sheets  were  subjected  to  a  pressure  of  280  psi.  at 
515°F  for  about  30  min.  and  then  cooled  rapidly  in  the  press. 
Before  the  test  specimens  were  formed  the  material  was 
given  a  second  anneal  by  heating  it  in  an  electric  oven  at 
275°P  for  15  min.  The  anneal  was  observed  to  cause  a  shrinkage 
of  5.7%   in  width  and  an  expansion  of  l.Sfo   in  length  measured 
at  the  surface  of  the  sheet. 

Specimens 

A  milling  cutter  was  used  to  cut  the  creep  specimens  from 
the  original  sheet .  For  polyethylene  and  Geon-4o4  the  flat 
side  of  the  specimen  was  perpendicular  to  the  plane  of  the  sheet, 
while  for  Kel-F  it  was  in  the  same  plane.  All  reduced  sections 
and  radii  at  the  ends  were  formed  on  a  shaper.  All  the  speci- 
mens were  finished  by  sanding  with  No.  0  emery  cloth.  The  di- 
mensions of  the  specimens  are  given  in  Fig.  1. 

Apparatus  and  Test  Procedure 

The  equipment  used  for  conducting  the  creep  tests  was  the 

(0) 
same  as  reported  previous ly^^''  for  creep  tests  of  five  laminated 

plastics . 

Creep  specimens  for  each  plastic  were  tested  simultan- 
eously under  a  constant  tension  load  at  several  different  values 
of  stress.   The  values  of  stress  chosen  were  evenly  spaced  from 
zero  to  a  maximum  value.  One  specimen  of  each  plastic  was 
tested  at  a  stress  of  zero  in  order  to  determine  the  magnitude 
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of  the  shrinkage  which  might  occur  due  to  gradual  changes  in 
moisture  content  or  other  time  dependent  phenomena,  such  as 
aging . 

All  creep  tests  were  conducted  at  a  constant  temperature 
of  77°  -  1°P  and  a  constant  relative  humidity  of  50±  2  percent 
for  the  entire  500  hours  of  the  test.   After  noting  the  initial 
extensometer  reading  Just  prior  to  application  of  the  load, 
readings  of  strain  and  time  were  taken  at  the  following  time 
intervals:  20  seconds,  then  1,  2,  3,  5,  7,  10,  12,  15,  l8,  30, 
k2   minutes,  then  1,  2,  3,  5,  7,  15,  l8  hours,  and  then  every 
2^  hours  to  about  500  hours . 

Further  information  about  the  apparatus  and  test  pro- 
cedure is  given  in  reference  (9). 

Results 

The  creep  curves  of  polyethylene  are  shown  in  Pig.  1. 
For  each  of  the  stresses  o  -   75,  150,  225  and  300  psi,  two  test 
specimens  were  used.   In  order  to  make  correction  for  changes  in 
relative  humidity  and  temperature,  an  additional  specimen  at 
zero  stress  was  used  as  a  control. 

The  creep  curves  of  Kel-F  (crystalline)  and  Geon-^04 
(annealed)  were  obtained  in  a  similar  manner,  i.e.  using  four 
different  stresses  and  a  total  of  nine  test  specimens  for  each 
of  the  plastics. 
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Evaluation  of  Creep  Constants 

In  order  to  describe  the  experimental  creep  curves 

by  Eq.  6,  it  was  necessary  to  evaluate  the  constants  n,  e', 

0  ,  m'  and  a„   for  each  of  the  plastics  by  using  the  correspond- 
o   o      m  «     ^ 

ing  experimental  creep  curves. 

The  following  procedure  was  employed  for  all  the  creep 
curves  of  a  given  material:  Rearranging  Eq.  4  and  taking 
logarithms  of  both  sides: 

log(e  -e^)  =  log  m  +  n  log  t  (7) 

Eq.  7  represents  a  straight  line  with  log  (e-e  )  as 
ordinate  and  log  t  as  abscissa.  Thus  a  particular  value  of 
€  was  selected,  such  that  when  log  (e-e  )  was  plotted  against 
log  t,  the  creep  data  fell  on  a  straight  line.  The  ordinate  of 
this  straight  line  at  t  =  1  hour  provided  the  value  of  m  while 
the  slope  of  the  straight  line  yielded  the  value  of  n. 

Values  of  €'  and  m'  were  obtained  by  the  use  of  Eq.  5. 

Appropriate  values  of  o^   and  a_  were  chosen  so  that  e^  vs 
^     ^  em  o 

sinh  — •  and  m  vs.  sinh  —  were  straight  lines.  The  slopes  of 

€  m 

these  lines,  then,  yielded  the  values  of  constants  €'  and  m', 

respectively. 

The  procedure  for  evaluation  of  these  constants  is 

illustrated  for  polyethylene .   If  strains  e  and  time  t  shown 

in  Fig.  1  are  plotted  on  a  log-log  graph,  the  data  fall  on  a 

curve  which  is  convex  downwards.   In  order  to  straighten  this 

curve  the  proper  choice  of  e^  was  made  such  that  when  log  (€-Gq) 

was  plotted  against  log  t,  the  data  were  on  a  straight  line. 
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The  data  for  each  different  stress  were  treated  in  this  manner 

and  are  shovm  as  straight  lines  in  Fig.  2.  The  slopes  of 

these  lines  are  the  same  and  yield  the  value  of  n,  while  the 

value  of  m  was  read  at  t  =  1  hr.   It  was  observed  that  both 

m  and  e^   increased  with  increase  of  stress  a.     Similar  diagrams 
o 

for  Kel~F  (crystalline)  and  Geon-koh   (annealed)  are  given  in 

Pig.  3  and  Fig.  4.  The  fact  that  the  data  are  nearly  on  a 

straight  line  indicates  agreement  with  the  data. 

For  evaluating  constants  €',  a  ,  m'  and  cr^,  of  poly- 

ethylene,  reference  may  be  made  to  Fig.  5.  Appropriate  values 

of  a^   and  a_  were  chosen  so  that  €^  vs  sinh  -— -  and  m  vs  sinh  -— 
€      m  o        a  0 

yielded  straight  lines.   In  order  to  increase  the  accuracy  of 
these  constants,  proper  scales  for  the  two  axes  were  chosen 
so  that  the  straight  lines  made  approximately  a  45°  inclination 
with  the  axes.  The  slopes  of  these  lines  yielded  the  values  of 
constants  e'  and  m'.  Similar  diagrams  for  Kel-F  (crystalline) 
and  Geon-4o4  (annealed)  are  given  in  Fig.  6  and  Fig.  J.     The 
values  of  the  constants  of  Eq.  6  are  given  in  Table  I  for  the 
three  plastics.   It  was  observed  that  for  Kel-F  (crystalline) 
the  value  of  n  at  the  highest  stress  a   =  56OO  psi.  was  much 
greater  than  for  the  lower  stresses.  A  much  higher  total  creep 
and  rate  of  creep  was  also  observed  at  this  stress.  This  fact 
is  shown  in  Fig.  5  where  the  slope  n  of  the  straight-line 
corresponding  to  stress  or  =  5600  psi  is  much  more  than  the 
constant  slopes  of  the  other  lines  at  lower  stresses. 

At  this  highest  stress,  a  noticeable  amount  of  crazing 
occurred  in  Kel-F  (crystalline)  which  probably  changed  the 
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structure  of  the  material  and  may  have  produced  the  increase 

in  creep.  The  stress  at  which  crazing  started  was  determined 

from  a  supplementary  creep  test  in  which  the  stress  was 

increased  at  intervals  from  2700  to  3000  to  3300  to  360O  psi. 

Crazing  started  after  19  hdurs  at  3300  psi. 

The  values  of  strain  €  computed  from  Eq.  6  are  plotted 

in  Fig.  1  for  Polyethylene.   These  theoretical  data  show  a  good 
comparison  with  the  observed  data  for  all  the  constant  stresses 
employed.  At  the  highest  stress  of  a  =  3OO  psi  the  theoretical 
values  seem  to  be  slightly  higher  than  the  observed  ones  for 
the  interval  of  time  t  =  0  to  t  =  250  hours .  This  discrepancy 
may  be  due  to  a  certain  amount  of  compromise  required  in 
selecting  the  constants,  see  Fig.  2  and  5^  or  to  errors  in 
observation  of  the  experimental  creep  curves . 

The  theoretical  values  of  strain  e  from  Eq.  6  were  com- 
puted also  for  Kel-F  (crystalline)  and  Geon-4o4  (annealed). 
For  these  plastics,  also,  an  equally  good  comparison  with  the 
observed  data  was  noticed  except  for  the  highest  stress 
a   =  3600  psi  for  Kel-F  (crystalline). 

20 -Second  Elastic  Strains 

The  extent  of  strain  in  a  plastic  after  it  has  been  under  a 
constant  load  for  20  seconds  is  designated  as  20-second 
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elastic  strain.  The  word  "elastic"  is  justified,  because 
the  contribution  of  the  time -dependent  strain  to  the  total 
strain  is  negligible  for  such  a  short  interval  of  time. 
The  20-second  elastic  strains  of  the  plastics  were  found  for 
each  of  the  several  constant  stresses  by  plotting  the  exper- 
imental strain  e  against  time  t  on  log-log  paper  up  to  a  time 
of  one  hour.  The  value  of  the  20-second  strain  was  then 
read  for  a  time  of  20- seconds  from  the  straight  line  which 
resulted  from  this  plotting. 

The  total  strain  (elastic  plus  plastic)  at  5OO  hours 
for  a  given  stress  o  can  be  read  from  the  creep  curves  such 
as  Pig.  1.  The  elastic  strain  at  20  seconds  and  the  total 
strain  at  5OO  hours  for  different  stresses  are  plotted  in 
Fig.  8  for  all  three  of  the  plastics.  It  was  observed  that 
the  appearance  of  these  curves  was  similar  to  stress-strain 
curves  in  tension. 

The  tangent  moduli  of  elasticity  obtained  by  measur- 
ing the  slopes  of  these  curves  at  about  zero  stress  is 
shown  in  Table  II.   A  comparison  of  the  moduli  of  elasticity 
obtained  from  the  20-second  elastic  strains  of  the  creep  curves, 
with  those  obtained  by  a  static  tension  test  at  the  same  temper- 
ature and  relative  hiimidity  would  be  desirable. 
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However,  the  values  of  the  moduli  of  elasticity  obtained 

from  the  20-second  elastic  strain  can  be  compared  with  those 

computed  by  a  stress-strain  relation  for  static  tension  from  the 

creep  Eq.  6.  This  relation  was  presented  in  a  previous  paper ^  ' 
in  the  following  form: 

n 


J. —11 
e  =  €•  slnh  ^  +p2^  slnh  f-]  .  m-  slnh  f-         (8) 


The  constants  €',  a^,  m»,  a^   and  n  in  Eq.  8  are  the  same  as 

o   e       m  ^ 

in  the  creep  Eq.  6,  while  V  expresses  the  strain  rate  of  the 

tension  test  curve. 

The  tangent  modulus  of  elasticity,  E,  was  obtained 

from  Eq .  8  by  taking  the  derivative  with  respect  to  stress,  as 

follows: 

n     1 

At  stress  a  =  0,  Eq.  9  yielded  the  relation 

E  =  ^  (10) 

o 

The  evaluation  of  E  from  Eq.  10  involved  the  determination 

of  the  constants  a^   and  €'.  These  constants  are  given  in 

e      o  , 

Table  I.  The  values  of  the  tangent  modulus  E  evaluated 
theoretically  from  Eq.  10  together  with  those  obtained  by 
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measuring  the  slopes  of  the  curves  for  the  20-second  elastic 
strain  in  Fig.  8  are  given  in  Table  II.   A  fairly  good  com- 
parison between  the  two  sets  of  values  was  observed. 

It  was  observed  that  the  value  computed  from  Eq .  10 
was  higher  for  each  material  than  that  measured  from  the 
slopes  of  the  curves  in  Fig.  8.  This  is  in  accord  with  the 
predictions  of  Eq.  9  as  described  in  reference  (lO)   In  this 
reference  it  was  shown  that  the  tangent  modulus  predicted  by 
Eq.  9  decreased  rapidly  with  increase  in  stress  for  stresses 
near  zero.  Since  the  slopes  measured  from  Fig.  8  were  inter- 
polated from  stress  values  substantially  above  the  origin 
the  values  of  the  moduli  should  have  been  less  than  predicted 
by  Eq.  10,  as  observed. 

Creepocity 
Equation  (4-)  for  creep  can  be  used  to  evaluate  the 
percentage  increase  in  strain  for  a  given  time  interval  (called 
creepocity).  See  reference  (9).  The  percentage  increase  in 
strain,  <t>,  from  t  =  0  to  t  =  t,  is  given  by 


(e^  +  mt^)  -  € 


1  ^      ~o 


2x  100  =f-(lOOt^)  =2_^^       (11) 


where 


^^      =  lOOt^  (12) 

(9) 

when  m  =  Gq. 'J't  is  independent  of  stress    ,  But  when  ; 

m  =  m'  sinh  ~-  and  e^  =  e^   sinh  |-  ,  4>-j^  is  a  function  of 
m  e 

stress  a   as  well  as  time  t^,  while  ^,  is  independent  of  stress 
and  dependent  on  time  t^  alone.   <^^  was  found  to  be  nearly 
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Independent  of  stress  for  five  plastic  laminates^^^ .  The 
"creepocity"  of  a  plastic,  when  independent  of  stress,  serves 
as  one  measure  of  its  serviceability  under  conditions  of  creep. 
For  the  materials  tested  in  this  investigation,  however,  the 
creepocity  4),  was  a  function  of  stress  since  e^  did  not 
equal  m.    ; 

If  the  percentage  increase  in  strain  is  determined 
from  t,=  20  sec.  to  tp  =  500  hrs.,  4"  is  expressed  as: 


0  =  [ 2 X  100  \   (t^  -  t^)  (15) 

I  1+  I  t^       J 
o  ~ 

The  values  of  <t)  determined  from  the  experimental  data 
and  given  in  Table  II  appear  to  show  no  systematic  increase 
with  increase  of  stress  for  a  given  plastic.  The  reasons' for 
this  are  probably  that  <J>  is  not  very  sensitive  to 
changes  in  stress,  and  insuiffic lent  accuracy  in 
determining  the  strains  experimentally.  However, 

it  was  noticed  from  Table  II  that  the  average  values  of  <t> 
obtained  theoretically  by  use  of  Eq.  15  compare  well  with  the 
average  value  obtained  from  the  experimental  data,  except  for 
polyethylene. 

It  was  observed  that  the  average  value  of  <J>  was  highest 
for  polyethylene  and  lowest  for  Kel-F  (crystalline)  which  means 
that  the  time-sensitivity  of  polyethylene  was  the  greatest 
while  that  for  Kel-P  was  the  least. 
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Rate  of  Creep 

The  rate  of  creep  V,  at  a  specified  time,  was  obtained 
from  Eq.  4  or  Eq.  6  by  taking  the  derivative  of  strain  €  with 
respect  to  time  t: 

V  =  ll  =  mnt""'^  =  (m-nt'^'^lsinh  §-  (l6) 

m 

The  rate  of  creep  at  500  hours  was  evaluated  by  using  Eq.  l6 
and  the  constants  m',  a  and  n  which  are  given  in  Table  I. 
The  results  are  shown  in  Table  II  and  are  plotted  with  respect 
to  stress  a  for  each  of  the  plastics  in  Fig.  9. 

For  a  given  stress,  the  highest  rate  of  creep  was 
observed  for  polyethylene  and  lowest  for  Kel-F  (crystalline). 
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TABLE  I.   CONSTANTS  OF  CREEP  EQUATION 
FOR  THREE  PLASTICS 


MATERIAL 

n 

'0' 
percent 

percent 

psi 

psi 

Polyethylene 

0.0890 

1.550 

0.597 

400 

185 

(Crystenme)  °-°^'^^*       °-^^°  °-°99    ^600     1475 


?lnn;a?ecl)    0-51°9    1.150     0.018    60OO     2100 


*  For  stresses  o  =  900,  I8OO,  2700  psi,  n  =  O.O872  while 
for  o  =  5600  psi,  n  =  O.I85. 
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FIG    2.        EVALUATION     OF  CONSTANTS     t„  ,  m  ,    AND  n     FOR     POLYETHYLENE. 
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AND  n  FOR    KEL-F  (CRYSTALLINE  ) 


a  '  4000 psi 


-  eo'0£22:  m'0.072 

•  -e„=0.820.m=0.0679 

a'iOOOpsi 

a  -e„'0.600:i7fO.0355 

cr'20O0  psi 

°-io'Ojaa:m'0.o/99 

•  -ecrOJ73;m'O.OI99 
a'lOOOpsi 

o-eff'0.184:  m'0.00785 
»-Sa'ai55:  m'0.00760 

n' 0Jt09    for   all  (r 
(m ,  Cq  in  percent} 


1.0 


10 


100 


1000 


Time,   hr. 
FIG    4       EVALUATION    OF    CONSTANTS     £„   ,  m  ,  AND  n    FOR    GE0N-404(ANNEALED), 
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EFFECT  OF  CRYSTALLINITY  AND  CRAZING,  AGING,  AND  RESIDUAL  STRESS 
ON  CREEP  OP  KEL-F,CAiWAS  LAMINATE,  AND  GEON  RESPECTIVELY 

1 
by  V/illiam  N.  Findley 

SMMARY 

Results  of  tension  creep  tests  at  77°  F  and  50  percent 
relative  humidity  are  presented  which  show  the  effect  of  the 
following  variables : 

(a)  the  degree  of  crystallization  of  rnonochlorotri- 
f luorethylene  (Kel-P) 

(b)  residual  stresses  remaining  after  fabrication  of 
sheets  of  polyvinylchloride  (Geon  4o^)  and 

(c)  aging  of  grade -C  canvas  laminate  at  77  F  and 
50  percent  relative  humidity. 

The  Kel-F  xms  tested  in  the  amorphous,  fully  crystalline 
and  intermediate  (commerical)  state.   The  effect  of  residual 
stresses  x-jas  shown  by  tests  of  samples  of  Geon  kok   before  and 
after  an  anneal  at  275°  F  which  resulted  in  shrinkage  of  the 
material.  The  canvas  laminate  was  tested  in  creep  before  and 
after  aging  at  constant  temperature  and  humidity  for  about  ten 
years . 
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CONCLUSIONS 

I,     Creep  of  Kel-P  was  markedly  influenced  by  crystallinity. 
The  most  crystalline  material  was  the  most  resistant  to 
creep . 

ft.  Crazing  of  Kel-F  was  found  to  result  from  creep.  At  the 
same  stress  the  amorphous  sample  showed  some  crazing, 
intermediate  samples  showed  none,  and  the  crystalline 
sample  crazed. 

3.  Creep  produced  crazing  in  the  crystalline  Kel-P  in 

19  hr ,  at  2210  psi,  and  at  56OO  psi  considerable  crazing 
was  present. 

4.  When  crazing  occurred  the  creep  strains  were  much  larger 
than  without  crazing, 

5.  Ilesidual  tensile  stresses  in  Geon  4o4  resulted  in 
greater  creep. 

6.  Aging  of  grade-C  canvas  laminate  produced  a  substantial 
reduction  in  creep. 
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INTRODUCTION 

The  effect  of  different  degrees  of  crystallinity,  residual 
stress  and  aging  on  creep  of  materials  has  not  received  wide 
attention.  The  only  work  on  these  topics  for  plastics  which 
has  been  reported,  as  far  as  is  known,  are  a  study  of  the 
effect  of  aging  on  creep  of  cellulose  acetate  (l)  and  a 
creep  test  method  for  determining  the  cracking  sensitivity 
of  polythelene  polymers  (2). 

Of  course,  most  engineering  materials  are  not  available 
with  a  wide  range  in  degree  of  crystallinity.  Metals,  for 
example,  are  largely  crystalline  and  contain  only  minute  but 
important  regions  of  amorphous  material  at  the  grain  boundaries 
Changes  in  structure  V7ith  time  are  found  in  other  engineering 
materials  besides  plastics  at  certain  temperatures,  such  as 
steels  at  high  temperatures.  Residual  stresses  may  be  found 
in  many  applications  of  many  materials.  Usually  residual 
stresses  are  relatively  low  in  magnitude  at  temperatures 
where  creep  becomes  significant  in  many  engineering  materials 
such  as  high  temperature  steels . 

The  effect  of  residual  stresses  in  plastics  has  often 
been  observed.  Means  of  reducing  these  stresses  and  the 
theory  and  measurement  of  residual  stresses  for  plastics 
have  been  considered  by  some  investigators.  Crazing  of 


Numbers  in  parenthesis  refer  to  the  list  of  references 
appearing  at  the  end  of  this  paper. 
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plastics  has  been  the  subject  of  several  investigations. 
For  linear  polymers,  crazing  has  been  described  as  a  mechanical 
separation  of  polymer  chains  or  groups  of  chains  under  the 
action  of  tensile  stress,  and  starts  in  a  region  where  the 
polymer  chains  are  oriented  at  right  angles  to  the  tensile 
stress .  References  to  these  topics  are  given  in  a  feature 
article  in  Applied  Mechanics  Revievjs  (5). 

MATERIALS 

Kel  -•  P  (Tponpchlorotrifluoroethylene)  of  grade  500  was 
obtained  in  sheets  of  I/8  in.  thickness  from  the  M.  W.  Kellogg 
Co.  This  grade  lAras  unplasticized  and  represents  the  highest 
molecular  weight  in  the  form  of  molding  powder.   The  high 
density  molding  powder  had  a  compression  ratio  of  2  to  1. 
The  sheets  were  compression  molded  under  a  pressure  of  2500  to 
3000  psi.  and  at  a  temperat^ore  of  480  to  500°  F.   Immediately 
after  pressing  a  sheet  was  quenched  in  v-iater  to  produce  the 
amorphous  sample.  Another  sample  v:as  cooled  between  platens 
which  had  cold  water  circulating  through  them  to  produce  the 
material  with  intermediate  crystallinity  (commercial  practice). 
The  crystalline  form  was  obtained  after  pressing  by  dropping 
the  temperature  quickly  by  means  of  steam  to  450°  P,  then 
cooling  to  room  temperature  in  the  press.   It  required  6  to  8 
hours  to  cool. 
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The  grade -C  canvas  lai^lnate  employed  in  these  tests  was 
the  same  from  which  test  specimens  were  obtained  for  previous 
tests  (4).  A  detailed  description  of  the  material  is  given  in 
reference  (4). 

Geon  4o4  (polyvinyl  chloride)  was  obtained  as  a  sheet 
9/16  inches  thick  from  the  Nixon  Nitration  Works .  The 
formulation  was  designated  VX-IOS  Natural  and  contained  no 
plasticizer,  3-3  percent  dibutal-tin  dlborate  as  a  stabilizer 
and  0.5  percent  stearic  acid  as  a  lubricant.  The  ingredients 
were  mixed  and  blended  at  280^  F  and  calendered  into  sheets 
0.02  in.  thick  on  rolls  having  a  temperature  of  300^  F, 
About  28  of  these  sheets  were  laminated  together  in  a  hydraulic 
press  to  form  the  final  sheet  and  remove  most  of  the  residual 
strain  resulting  from  the  calendering.  The  sheets  were  sub- 
jected to  a  pressure  of  280  psi  at  315°  F  for  about  30  min. 
and  then  cooled  rapidly  in  the  press. 

SPECIMENS  AND  TEST  FROCSPmE 

The  specimens  and  procedure  employed  for  the  creep  tests 
Vfere  the  same  as  previously  described  (4)  except  that  the 
specimens  for  creep  tests  of  Kel  -  P  and  Geon  4o4  had  4-inch 
gage  length  and  the  specimens  of  Kel  -  F  were  the  thickness  of 
the  sheet  instead  of  5/32  in.  thick. 
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RESULTS  AND  INTERPRETATION 

Figure  1  shows  the  results  of  the  creep  tests  of  Kel  -  P 
in  three  different  degrees  of  crystallization.   It  was 
observed  that  the  greatest  creep  occurred  in  the  amorphous 
and  the  least  in  the  crystalline.  Both  the  time -independent 
(elastic)  strains  and  creep  strains  were  much  larger  in  the 
amorphous.   The  total  strain  at  a  stress  of  2700  psi  after 
500  hours  under  a  constant  load  was  5.3^  2.8,  and  1.5  percent 
for  the  amorphous, intermediate,  and  crystalline  respectively. 
Recovery  following  creep  of  the  amorphous  and  intermediate 
samples  is  also  shown  in  Pig.  1. 

Creep  of  the  crystalline  material  is  shown  for  four 
different  stresses  in  Fig.  1.   It  was  observed  that  creep 
at  the  stress  of  3600  psi  was  much  greater  than  would  be 
expected  from  the  creep  of  the  crystalline  material  at  lower 
stresses. 

Examination  of  the  specimens  disclosed  that  crazing 
occurred  in  the  crystalline  material  at  the  stress  of  36OO 
psi.   It  was  not  known  whether  crazing  occurred  in  the 
lower  stresses.   Crazing  was  not  observed  in  the  material 
of  intermediate  crystallinity,  but  slight  crazing  was  ob- 
served on  one  side  of  the  amorphous  sample. 

This  suggests  that  the  larger  creep  in  the  sample 
tested  at  360O  psi  was  probably  due  to  the  combination  of 
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crazing  plus  creep.   In  order  to  determine  at  what  stress 
crazing  started  in  the  crystalline  material,  a  creep  test 
was  performed  in  which  the  stress  was  increased  at  intervals 
from  1610  to  1910  tc  2210  to  2510  psi.   It  was  found  that 
no  crazing  occurred  until  about  19  hours  after  the  load  was 
increased  to  2210  psi.  At  this  time  the  strain  increased  from 
1.06  to  1.27  per  cent.  From  this  time  on,  a  count  was  taken 
of  the  number  of  transverse  crazing  marks  per  inch.   In  one 
inch  the  number  of  marks  increased  in  succeeding  24  hour 
periods  from  7  to  15  to  125  to  I30. 

Results  of  two  creep  tests  of  Geon  4o4  at  40CC  psi  are 
shown  in  Fig.  2,  together  vjith  recovery  follov7ing  removal  of 
load.   One  of  these  specimens  was  prepared  from  the  sheet  of 
material  as  received  from  the  molder.  Tne  other  s^recimen  was 
prepared  from  the  same  sheet  v.^ith  the  same  orientation,  after 
the  sheet  of  material  had  been  annealed  to  remove  residual 
stresses.   The  material  was  annealed  by  placing  it  on  a  glass 
plate  and  heating  in  an  electric  oven  at  275°  ^  for  15  minutes 
and  cooling  slowly.  Measurement  of  reference  marks  on  the 
surface  before  and  after  the  annealing  indicated  a  shrinkage 
of  5.1  and  0.8  percent  lengthvjise  and  crossv;ise  respectively 
of  the  original  sheet  (and  direction  of   the  calendaring). 
Greater  shrinkage  at  the  center  of  the  sheet  than  at  the  sur- 
face was  indicated  by  the  concave  appearance  of  the  edges  of 
the  sheet. 
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As  shown  in  Fig.  2  the  as -received  sample  showed  more 
creep  than  the  annealed  sample,  about  1.48  and  1.59  percent 
creep  respectively.  Apparently  the  tensile  residual  stresses 
locked  in  the  sample  made  the  material  more  susceptible  to 
the  action  producing  creep. 

Creep  tests  at  two  different  stresses  are  shown  in  Pig. 3 
for  a  grade-C  canvas  laminate.  The  two  tests  shown  at  each 
stress  were  tested  with  an  elapse  of  tim.e  of  about  10  years 
between  tests.  When  tested  in  August  19^3  the  laminate  was 
about  one  month  old.  It  has  been  stored  at  a  constant  tem- 
perature of  77°  F  and  50  percent  relative  humidity  continu- 
ously since  that  time  except  during  machining  of  specimens. 

The  results  shown  in  Fig.  3  indicate  that  the  creep 
resistance  has  improved  with  age.  The  creep  at  500  hours 
was  reduced  with  age  from  1.17  to  0.99  percent  at  5700  psi 
and  0.^2  to  0.57  percent  at  2850  psi. 
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